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Preface 


This monograph deals with the problems of mathematical physics which 
are improperly posed in the sense of Hadamard. 

The first part covers various approaches to the formulation of improperly 
posed problems. These approaches are illustrated by the example of the classical 
improperly posed Cauchy problem for the Laplace equation. 

The second part deals with a number of problems of analytic continuations 
of analytic and harmonic functions. 

The third part is concerned with the investigation of the so-called inverse 
problems for differential equations in which it is required to determine a dif- 
ferential equation from a certain family of its solutions. 


Novosibirsk June, 1967 M. M. LAVRENTIEV 
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Chapter | 


Formulation of some Improperly Posed Problems 
of Mathematical Physics 


§ 1. Improperly Posed Problems in Metric Spaces 


The notion of correctness* introduced at the beginning of our century by 
the French mathematician HADAMARD plays an important role in the investiga- 
tion of the problems of mathematical physics. One often says that a problem is 
solved if its correctness is established. 

Various authors present notions of correctness which coincide in their 
essence but differ in details. We give one of the possible definitions of correct- 
ness which is convenient for our aims. 

Let ®, F be some complete metric spaces, and let Ag be a function with the 
domain of definition ® and the range of values F. Consider the equation 


Ao=f (1.1) 

Let us point out that most problems of mathematical physics can be reduced 
to the investigation of the solution of equation (1.1) with a given function A 
and right-hand side f. We say that the problem of solving (1.1) is properly 
posed if the following conditions are satisfied: 

1) The solution of (1.1) exists for any fe F. 

2) The solution of (1.1) is unique in ®. 

3) The solution of (1.1) depends continuously on the right-hand side f- 


In other words, the problem of solving (1.1) is properly posed if there exists 
a function Bf defined and continuous over all of F, which is inverse to the 
function Ag. 

Linear problems are most often considered in mathematical physics. In 
this case ®, F are BANACH spaces, and 4 is a linear operator. The BANACH 
spaces ®, F in concrete problems are the known functional spaces C', L,, Ww, 
A p S,,--. with the carriers in some n-dimensional space of the independent 
variables or on any part of the spaces of independent variables. 


* Translator’s note: Hereafter we shall refer to the correctness or incorrectness of 
problems which are respectively properly posed or improperly posed. 


1 Springer Tracts, Vol. 11, Lavrentiev 
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The first requirement of correctness is that the problem should not be 
overdetermined, and superfluous conditions should not be imposed. 

The second requirement is that the solution be unique. 

The third requirement of correctness, continuity of the inverse function Bf, 
arises from the fact that in the real problems of mathematical physics the right- 
hand side of equation (1.1) is obtained from measurements made with the aid 
of actual instruments and is therefore known only approximately. Therefore, 
it has been felt for a long time that if at any point f the function Bf is dis- 
continuous, then the solution g cannot be uniquely recovered from the right- 
hand side f. 

HADAMARD introduced the notion of correctness by giving an example of 
an improperly posed problem which became classical and was included in 
most text-books on mathematical physics. The example is the CAucHY problem 
for the LAPLACE equation. It is well-known that the solution of this linear 
problem does not depend continuously on the data obtained from any of the 
functional spaces mentioned above. On the basis of this HADAMARD concluded 
that CaucnHy’s problem for the LAPLACE equation and, in general, all problems 
exhibiting a similar dependence of the solution on the right-hand side, do not 
correspond to any real formulations, i.e., they are not problems of mathematical 
physics. 

It was discovered later that HADAMARD’s conclusion was erroneous and 
many real problems of mathematical physics lead to problems which are 
improperly posed in the sense of HADAMARD. In particular a number of im- 
portant problems of geophysics lead to the CAucHy problem for the LAPLACE 
equation. 

There are still further examples of linear as well as nonlinear improperly 
posed problems which are important in the applications, namely the solution 
of the heat equation for negative time and CAucuy data on the boundary, the 
nonhyperbolic Caucuy problem for the wave equation, inverse problems of 
potential, and a number of inverse problems for differential equations. 

At present there exists a number of approaches to the investigation of 
improperly posed problems. We will now explain them using the above classical 
Caucny problem for the LAPLACE equation as an illustration. 

We consider one of the simplest versions of the CAucuy problem for the 
LAPLACE equation. Let u (x, y) be a twice continuously differentiable function 
of the variables x, y in the rectangle 


O<x<zn, O<y<H (D) 
satisfying the following conditions 
Au=0, (x,y)eD (1.2) 


0 
apt D=vO.y) =u, y)=0- 
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It is required to determine its values on the segment y = h, 0 < x < a by 
its values on the segment y= 0,O< x < x. 


The problem stated is equivalent to the solution of equation (1.1) where 
y, f stand for the functions of the variable x 


g=u(x,h) 
f =u(x,0). 


A is a linear integral operator with the kernel 


>Y cosh 'kh-sinkx:-sinké. 


2 
Tl y=1 


K, (x, f= 


The functional spaces ®, F considered in the given example are the HILBERT 
functional spaces with the scalar products of the form 


(91, 92)= » Ae) @4 p;e® 


0 (1.3) 
i Sd= Yi mst fr S,;€F 


where Cy and f;* are the FOURIER coefficient of the functions g; and f;, and A, 
and 4, > 0 are some sequences. 

We note that spaces of such form are, in particular, the spaces wi with the 
norm 
ond 
Ox4 ‘ 


I 
lvlly = > C; 
J=0 


L2 


It may be easily seen that the solution of (1.1) in the case under considera- 
tion is unique, but for its existence and continuity it is insufficient to let ® and 
F be some W, with finite /. 

To make the problem properly posed one can introduce a sufficiently 
‘‘strong’’ norm in the data space F or a sufficiently ‘‘weak”’ one in the space ®. 
It is quite evident that the solution of (1.1) is correct in the sense of our defini- 
tion, if the sequences A,, 4, in (1.3) satisfy the inequality 


Ae<C by, (1.4) 


where C is a constant. 


The inequalities (1.4), in particular, are satisfied by the following sequences 
Ay, 4, Corresponding to the above possibilities 


A=e™, by=1 


5 
A=1, mae’. ae 
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However, the above approaches have the following defects. In the first 
case the statement does not cover the range of problems in which the errors 
in the data may be considered to be small only in the “‘ordinary”’ functional 
spaces; in the second case solutions of the problem may be objects which are 
not functions in the usual sense and do not correspond to physical reality. 

We now formulate some approaches to the question of correctness of 
problems of the type under consideration which are free of the above defects 
and, in our opinion, quite natural from the standpoint of applications. 

The first approach consists of changing the notion of correctness, namely, 
to one having requirements different from 1), 2), and 3). 

We now state the requirements. In addition to the spaces ® and F and the 
operator A, let there be given some closed set Mc ®. 

We call the problem for the solution of (1.1) properly posed according to 
TykuHoNov if the following conditions are fulfilled. 


1) It is a-priori known that the solution @ exists for some class of data and 
belongs to the given set M, pe M. 
2) The solution is unique in a class of functions belonging to M. 


3) Arbitrarily small changes of the right-hand side of f which do not carry 
the solution @ out of M correspond to arbitrarily small changes in the solu- 
tion @. 

We denote by MZ, the image of M after the application to the space @ of the 
operator A. 

Requirement 3) can be restated in the following manner, 

3) The solution of equation (1.1) depends continuously on the right-hand 
side fon the set M4. 

If M is a compact set the following statement holds (see [17]). 

If equation (1.1) satisfies the requirements 1), 2) of correctness due to 
TYKHONOV, then there exists a function « (t) such that 

a) « (t) is a Continuous nondecreasing function with « (0) = 0. 


b) for any @;,@2 € M satisfying the inequality 


@(Ag,,Ag2)<E 
the following holds 


Q(91, P2)<a(e). 


Thus, the requirement of continuous dependence 3) is satisfied if 1) and 2) are 
satisfied. 


We note that, if a problem is properly posed according to TYCHONOV and 
we replace the metric spaces ®, F by their supspaces M, M 4 then the problem 
becomes properly posed in the usual sense. 
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The necessity of examining spaces ©, F together with M, M 4 is due to the 
fact that in real problems the errors committed in the determination of the 
right-hand side f usually lead to f outside of M,. The consideration of the 
problem according to TYKHONOV’s formulation gives the possibility of construct- 
ing an approximate solution with a certain guaranteed degree of accuracy in 
spite of the fact that an exact solution of (1.1) with approximate data either 
does not exist at all or may strongly deviatc from the “‘truc”’ solution. 

In the Caucny problem for the LAPLACE equation we consider as the sct M 
the sets defined in the following manner: we denote by A, an integral operator 
with the kernel 


2 pig 2 ; 
Ky, &, =) cosh”! kh,sinkx-sinké 
0) 
and let 


geM(I,h,,C) 


p=A,,Y, wily w<C, (h,=H—h). 


We now give evidence that our formulation is natural for the case in which 
the CAucny problem for the LAPLACE equation is used for the solution of the 
geophysical problem of interpreting the gravitational or magnetic anomalies. 
The problem of interpretation of geophysical data is as follows: the characteris- 
tics of some physical field interacting with the inner layers of the earth are 
measured on the earth’s surface. Certain characteristics of the structure of 
these layers are required to be determined. In interpreting the constant magnctic 
and gravitational fields there often arises the following situation. An upper 
layer of the earth’s crust consists of uniform sedimentary rocks. The thickness 
of the sedimentary rock layer may be estimated from below by some constant 
M. The anomalies observed on the earth surface, i.e., the CAUCHY data for a 
harmonic function, are generated by some bodies lying at a depth exceeding 
and the ultimate aim of an observer is to determine these bodics. This problem 
is rather complicated; its solution is not unique if one proceeds only from 
measurements of a field on the surface. However, by proper evaluation of the 
general geological situation and with well chosen supplementary hypotheses, 
the observer often succeeds in finding solutions with satisfactory accuracy. 

It is important in this problem to make good estimates on the number of 
the bodies, their approximate position, and their size. For the case in which 
the average depth of embedding of the bodies is substantially less than their 
horizontal dimensions and the spacing between them, the distinct extrema of 
the anomaly correspond to distinct bodies. However, if the average dcpth of 
imbedding roughly coincides with or excceds the characteristics horizontal 
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dimensions, then there is no such correspondence. It is possible for two bodies 
to cause one extremum situated in the gap between them. In order to obtain 
the case in which the depth of embedding is less than the horizontal dimensions 
one uses the solution of CAUCHY’s problem for the LAPLACE equation. The 
existence of the solution is ensured by the fact that the field to be measured is 
generated by real objects and by the fact that in the region of evaluation, the 
depth of which does not exceed H, singularities of the field do not occur. The 
fact that the solution belongs to the set M is guaranteed by the fact that in the 
earth’s crust one does not encounter bodies of excessive density or with magnetic 
intensity exceeding some completely determined constants. For example, when 
u(x, y) is a gravitational potential the above considerations guarantee tliat 
the solution belongs to a set M (2, A,, C), in the case when u is a component 
of the gravitational field stress—// (1, 4,, C), and in the case when u is a 
component of the magnetic field stress— M (0, /,, C). 

This approach to the question of the formulation of improperly posed 
problems was originally proposed by TyKHONOv [17]. A systematic investiga- 
tion of various problems, properly posed in the sense of TYKHONOV, was carried 
out in [28], [31], [50], [56]. A detailed description of the definitions as well as 
the investigation of a number of problems which are properly posed in the 
sense of TYCHONOV are contained in the monograph [35]. 

Another possible approach to the study of improperly posed problems is 
to let the notion of correctness remain unchanged, but change the notion of a 
“solution” of the problem. Besides the spaces ® and F and the operator A let 
there be given, as in the preceding case, some set M. By a quasisolution of 
equation (1.1) we mean an element @ defined by the relation 


0 (Ag, f)=min e(Ag, f). 
geM 

We note that the problem of finding a quasisolution is one of nonlinear 
programming. From the theorems of nonlinear programming one obtains the 
corresponding theorems on the correctness of the problem of finding a quasi- 
solution. 

The problems which may be regarded as properly posed according to 
TYKHONOV are those for which the notion of quasisolution may be suitably 
introduced. In some cases the quasisolution may be non-unique, but when the 
problem is properly posed due to TYKHONOV, any quasisolution constructed 
using an approximate right-hand side will be one of the best versions of an 
approximate solution of the problem. 

The notion of a quasisolution was introduced in [45]. Even though, as 
previously mentioned, the range of applicability of this notion coincides with 
that of the notion of correctness due to TYKHONOV, the notion differs by offering 
greater clarity and simplicity. 
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The most general approach to improperly posed problems was given in the 
recently published works of TyKHONOV. Following [18] we define a regularizor 
for Eq. (1.1) to be any one-parameter family of operators B, having domain of 
definition F and range of values © and satisfying the following conditions 

1) for any t > 0 the operator B, is defined and continuous over all of F. 

2) for anype®@: 

lim B,Ag=@. 
t70 

We call the problem of solving (1.1) regularizable if there exists a regularizor 
for (1.1). It may be easily scen that the regularizor permits us to construct a 
solution with a guaranteed degree of accuracy according to the approximate 
right-hand side. In fact, let f be the right-hand side of (1.1) with an error e, 
1.€., 


o(f, fe. 


We denote 
G.=Bf; 9,=BAg 
and estimate the quantity 0 (@,@,). 
By the triangle inequality 


2(9,6)<0(9,9,)+0(G. §,)- (1.6) 


It is clear that 


2(P, 9.) Sm (t) 


Q (9, @,) < X> (t, &) (1.7) 


where «, (t) is a function characterizing the rate of the convergence of B,Ag 
to, and «, (tT, €) is the modulus of continuity of the operator B, at the point 
Aq. The functions «, (t) and «, (t, €) for fixed t are continuous and monotonic, 
and a, (0) = «, (t, 0) = 0. 

Substituting (1.7) into (1.6) and setting t equal to a root of the equation 


Gt, (tT) =o (T, €) (1.8) 
where &, (7), is the inverse of «, (1), we obtain 
o(y, P)S2 Hy (t). (1.9) 


It is evident that the root of equation (1.8), and therefore, the right-hand 
side in the inequality (1.9) tends to zero as € > 0. 

In contrast to the first two approaches, in the notion of regularizability we 
consider no additional sets other than the basic spaces ®, F. This accounts for 
the greater generality of this approach. Apparently, all the problems of mathe- 
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matical physics connected with real phenomena are regularizable. (HADAMARD 
regarded all such problems to be properly posed.) In [22] the regularizability 
of a rather wide class of equations (1.1) has been proved. 

In conclusion we point out the relation between the notion of a quasisolu- 
tion and a regularizor. Suppose the solution of (1.1) is unique. Let M, be a 
one-parameter family of compact sets (t > 0) satisfying the requirements 


M,,¢M,,, ™|<t,U M,=®, 
t>O 
and let B, be a continuous operator with domain of definition F such that B, f 
is a quasisolution of (1.1) with the set M = M,, It is evident that the operator 
B, is a regularizor for cquation (1.1). 
Besides the above three approaches to improperly posed problems, an 
approach related to the theory of probability has been developed to some 
extent. The next section is devoted to this approach. 


§ 2. A Probability Approach to Improperly Posed Problems 


The theory of probability was first utilized in the investigation of improperly 
posed problems in a paper by KHALFIN and SUDAKOV [66]. The authors restricted 
their investigation to the classical CaucHy problem for the LAPLACE equation. 
The main content of the paper is as follows. 

The Caucuy data are assumed to consist of a sum of two items, the “exact 
data” and the “error”. The “error” is a realization of a stationary random 
process. The question to be investigated is, under what conditions on the 
autocorrelation function of the error will this error insignificantly affect (in the 
average) the solution of a problem. It turns out to be sufficient that the spectral 
function of the error should be small in the HILBERT space with the exponential 
weight. 

KHALFIN and SUDAKOV’s approach is analogous to a “‘determinate’’ approach 
in which the class of permissible initial data (see Section 1) is narrowed. We 
present here a statistical approach similar to the determinate one in which the 
notion of correctness according to TYKHONOV or of quasi-solutions is introduced. 
In this case we shall use the notions and definitions of information theory 
presented in the paper by DoBRUSHIN [49]. 

As it was noted in Section 1, most of the problems of mathematical physics 
reduce to the examination of equation 


Ag=f (1.1) 


where ¢, fare elements of the metric spaces ®, Fand A is acontinuous operator. 
In practice the right-hand side is, as a rule, a result of observations made with 
the aid of real physical instruments and therefore cannot be considered to be 
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given with absolute accuracy. However, in practice it is often the case that a 
certain statistical distribution of errors in the definition of the right-hand side 
of (1.1), and the statistical distribution of solutions @ are given. Motivated by 
this we propose the following approach. 

Let ©, F be measurable spaces (sec [49]), i.e., in the spaces ®, F besides the 
metrics, there are also given some o-algebras S, S, of subspaces of the spaces 
@®, F and the operator A is measurable. In addition let 2 be a probability 
space, 1.e., a measurable space with o-algebra Sy and with probability measure 
P { } given on this o-algebra. The clements ¢, f in (1.1) will be considered to 
be random variables with their values in the spaces ?, F. Morcover, we consider 
also a random variable f with values in F—the right-hand side of (1.1) with 
an error—and the random variable (f/, f ) with values in F x F, the product 
of F with itself. Let B be an operator with domain of definition in F and the 
range of values in ®, and let @ = Bf, where, @ are random variables in the 
space DP x @, 

The problem of constructing an approximate solution to (1.1) consists of 
constructing such an operator B that the distance between ¢, @ be sufficiently 
small where the notion of smallness should take into account both the metric 
of the space ® and the distribution of random variables (Q, @). 

As an estimate of the deviation of the approximate solution @ from the 
exact solution g one may take, for instance, the quantity mo(g,@)* or the 
quantity P {0 (g, @) > e}. Corresponding to these, we introduce two notions 
of optimal solution of equation (1.1). 


1. The function @ = Bf is called an optimal solution of (1.1) if 
mo(,@)=min. 
2. The function G = Bf is called an Sct € solution of (1.1) if 
P{ol9, 9) = eam. 


If the above estimates of deviation of @ from¢ do not exceed the number 6 
we call @ a solution of (1.1) with accuracy 6 in the sense 1 or 2, respectively. 
This statement is in its content close to the decoding problem of information 
theory (sec [49)). 

We note that since in our statement the operator A is determined, the 
distribution of the random variable g in © uniquely induces the distribution 
of the random variable f in F while the joint distribution of f, f in FX F 
induces the joint distribution of the random variables ¢, f, in the space ® x F. 
Since f in the statement is considered to be given, we may consider the dis- 


* mis the mathematical expectation, the random variable 0 (9, ¢) is assumed to be 
measurable. 
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tribution in ® induced by a joint distribution of g, f at fixed f. If the dis- 
tribution is to be given, the above formulated problem of constructing optimal 
solutions and solutions with prescribed accuracy become classical problems 
of the theory of probability. However, the practical construction of the joint 
distribution ofg, f may involve considerable difficulties, and so in most cases 
it is expedient to choose simpler although less accurate methods. 

The determinate approaches in which the notion of correctness according 
to TyKHOoNov or of a quasisolution are introduced, may be regarded as parti- 
cular cases of the proposed approach. In fact, let the set M7 introduced in the 
previous section be measurable, and 


P{peM}=1. (1.10) 


Clearly the equality (1.10) may be regarded as a priori information that 
the solution @ exists and belongs to the set M. 

Proceeding from the above ideas, we now give some results relating to 
optimal solutions and to solutions with a prescribed degree of accuracy, 
obtained under some additional assumptions. 

1. Let F, ® be HiILBerT spaces. We denote by pra random variable in ® 
whose distribution is induced by a joint distribution of g, f for fixed f. One 
can easily see that in this case the optimal solution in the sense of 1 is given 
by the formula, 

b= |erdo (1.11) 


In fact, for any 
m||6—pyll?=IG|l?-2(G, J wzpdeo)+ J |lp,ll? do 
2 2 


from which follows the validity of our assertion. 


2. Let the solution of (1.1) be unique, and let (1.10) hold, where M is some 
compact set and M, is the image (in the space F) of the set M. Then, as noted 
in the previous section, the operator A~! on a set M, is uniformly continuous. 
Its modulus of continuity will be denoted by « (t). 

Next let 


P{o(f, f)>d}=0 (1.12) 


where 6 is a constant. 
It can be easily seen that if ¢ in definition 2 of an optimal solution satisfies 
the inequality 
E<a(0) (1.13) 


then the quasisolution of (1.1) defined previously with the set M is an optimal 
solution in the sense of 2. 

In fact, let @ be a quasisolution of (1.1) with the righthand side f and the 
set M. 
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Then, by virtue of (1.10), (1.12), (1.13) and the definitions of quasisolution, 
and modulus of continuity « (z), we get 


P{o(9,6)>e}<P {o(f, f)>5}=0 
which means that the quasisolution @ is ¢ optimal in the sense of 2. 


3. Let B, be a regularizor of (1.1). 
We estimate the accuracy, in the sense of 2, of the solution of (1.1) with the 
aid of the regularizor B,. We assume, in addition, that the operator B, is uni- 


formly continuous on F for any 4 > 0, and denote by « (e, 2), «, (t), B (t) the 
functions 


a(e,A)=P {0 (B,AQ, y)>8} 


%,(t)= max o(B,u, B,v) 
e(u,v)<t 


B(t)=P (oF, f)>7} 


By virtue of the above assumptions and the known properties of the random 
variables 


P (0(B,F.0)<e}2max| PY o(B,F, 8, nsil: PloiBss, ost] (1.14) 


Py owss ,<shaPods)<ai(§)}=1-0] «(5)| (1.15) 


where o; (t) is the function inverse to a, (t). 
Substituting (1.15) in (1.13) we get 


P {0(B,f, 9) >e}<max LE (=) 5 2)] (1.16) 


We note that the estimate (1.16) makes possible an optimal choice of the 
parameter 2 of the regularizor B,. 


4. We now give the estimates of the effectiveness of a particular regularizor 
for the Caucny problem for the LAPLACE equation considered in the preceding 
section. Let F and @ be the space L, of functions defined on the segment [0, 7] 
and A, be an integral operator 


Avo=1Kr (x, )p(E) de 


with the kernel K,, equal to 
2 z ‘ 
K, (x, as cosh” ' kh-sinkx:-sin ké. 
0 


Next, let the distribution of the random variables g, f, f be such that 
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= 


P{\|A,, gl >t}=(pt? +1)" 


. (1.17) 
PUIS—fll>pH@r+h 
As a regularizor B, we consider an integral operator with the kernel 
9 . : 
0, (x, d)== ¥) coshkh(t + Ae)"  sinkx-sinké h,>h, 
TO 
Onc can show that the operator B, satisfies the following relations 
B,|| <a"? -r, (hh 
[Bl 1 (h, ha) (1.18) 


Ajo] 271 (hy, ha): (EB, A) | 
where 

(A, hy) = h"!"? (hy — hy" Ay! 

r> (hy, hy) =(hy—h,) 02 Pe Ae hy 


From (1.17) and (1.18) we get 


P{|B, fF —Nil>e}sP (IBF —fll>e}=(Qarn”- er 741) ¢~— 
P{|B,A—E)gll>e} <P {A,,\ ol] >a rye} = (P22 rt = 
from which it follows that 


P{|B,f —oll>2e}<max[(q?7" ery? +1) 1s (pea Pr $I)” ']. 
(1.19) 


It is evident that the first expression in the right-hand side of (1.19) mono- 
tonically decreases as A increases, while the second one monotonically increases. 
This means that in the sense of the estimate (1.19), the value of the parameter 
of the regularizor will be optimal when the two expressions coincide. By 
equating the expressions in the right-hand side of (1.19) we obtain 


i oe! PVs h2/(h thy) 
0 q & 


P {Bi Jf - ¢| >26} <(puerter. ghillathn), phalththy), rs 4 i! ; (1.20) 


From (1.20) it follows that for any 6, e > 0 there can be found a sufficiently 
large g— the index of the ‘“‘quality of the measurements’’— so that the estimated 
probability will be less than 6. 

We also note that the regularizor B, used by us also permits optimality 
with respect to the parameter /,, but the analytical determination of the optimal 
hy leads to extremely complex calculations. It seems that the parameter should 
be chosen on the interval [H, 2 7]. 


Chapter II 


Analytic Continuation 


The problem of analytic continuation is investigated in this chapter. Esti- 
mates of “‘stability”’ for the solution of some problemis are given, and concrete 
solutions are constructed. 


§ 1. Analytic Continuation of a Function of One Complex Variable from 
a Part of the Boundary of the Region of Regularity 


The problem was first considered by CARLEMAN in [3]. The results were 
obtained in [5]. 


I. The First Problem. 


Let f (z) be an analytic function, regular and bounded in some bounded 
domain D so that 
If(D|<M zeD (2.1) 


I’ is the boundary of D; J” isa part of J’and J” = J’ — I”. Let the values of 
f ( on I” be known, and suppose it is required to determine f (z) in some 
part of D. We now prove a theorem characterizing the stability of the solution 
of the problem. 


Theorem 1. Let f (z) on the curve /” satisfy the inequality 


If(z)lse (zer”) (2.2) 
Then the inequality 
If (z)|< MP7? +6? (2.3) 


will also hold where w (z) is the harmonic measure of the curve J” with respect 
to the point z and the domain D. 
For the proof we consider the function 


e(z)=In|f (z)]. 


It is known that @ (z) is a harmonic function. From (2.1), (2.2) it follows that 
the function ¢@ (z) satisfies the inequalities 
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y(z)<Ine zel" 
y(z)<InM zel” 

from which it follows that 
o(z)<w(z)Ine+[1—w(z)]InM (2.5) 


The inequality (2.3) to be proved follows from (2.5). 

The problem we are considering can be reduced to an integral equation 
of the first kind and can be solved by one of the general methods. In fact, we 
denote by D’ the part of D in which f (z) is to be determined. We denote by 
P, Q the ends of the curve J” and consider some points P;, Q; (j = 1, 2) lying 
on I” sufficiently close to the points P, Q, respectively. We denote by Lr; the 
part of J” lying between the points P;, Q;, and let the location of the points 
P;, Q; be such that I’; € Jj (see Figure !). 


(2.4) 


We draw a curve J’; through P,, Q, which lies inside D so that the domain 
D’ lies inside the domain D, bounded by the curves I I ; D’e Dj. It is evident 
that in order to determine the function f (z) for any zé D’ it is sufficient to 
determine f (z) on the curve J’, in view of the Caucny formula at ze D’ 


1 ff @) 
=—— ——d 
f@) Piette $ 
and the f (z) values given on I°{. 
We denote by f, (z) the limiting values on the curve J of the Caucny-like 


integral of the function f(z) on the are I” as z tends to Js from inside of D; 
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and by /f; (z) we denote the values on the curve J’; of the CAucHy-like integral 
of the function f(z) on the arc J” 


fa@mzalim (© 


Leo, 
zeD 


1 
dC; h@=sx5 [A Dat, zel,. 


The limiting values of the CaucHy-like integral are calculated by the familiar 
SOKHOTSKY formulas, so that for an investigation of this problem we may 
suppose that the functions f5, f, are given. 

We denote by 4 the operator transforming the arbitrary complex function 
@ (z) on I’ to the function w (z) given on Ij by the formula 


Wz asa J Peat eres 


It can be easily seen that the validity of the equation 


Ag(x)=f(z)-fa(z)  zel2 


is necessary and sufficient for the value of the function @ (z) to coincide with 
the values on J’, of the unknown function f(z). Thus, the problem of deter- 
mining the function f(z) on J’, is equivalent to that of solving the linear 
operator equation of the first kind 


Ag=f—fr. 
Then by the Caucny formula for any ze I, 
, 1 ze) 
$@=55 Set h@ 


from which it follows that the function f(z) on J’; is equal to 


f (z)=BW(z)+f, (2), 


where W (z) is a function defined on [”’ 
IW (ZI <M, 


(M is a constant depending on D), and B is the operator transforming the 
complex function W (z) given in I” to the function BW (z) given on I; by the 
formula 


1 ,W@, 
(2) 20; ,C-Z a4 ae 
It is clear that B is a linear completely continuous operator and therefore 
in the problem of solving the equation 
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Ag=f—-f2 
the conditions under which one treats linear operator equations of the first 
kind are met. 
We now consider another method of concrete solution based on the integral 
formulas of CARLEMAN type. We call a CARLEMAN function for the domain D 
and the curve J” any function G (z, ¢, 6) having the following properties 


1) G(2,,5) =p O(,6,8) 


where G (z, ¢, 6) is an analytic function of the variable ¢, regular and bounded 
in the domain D. 


2) the function G (z, ¢, 6) satisfies the inequality 


f IG(z,G dy |dt|<6. 


yt 
We will construct the concrete solution for the problem with the aid of a 
CARLEMAN function. Let the values of f(z) on the curve J” be known with 
accuracy ¢, i.e., a function f, (z) is known on the curve J” such that 


lf.(z)-f(zZlse, zer’. 
We denote by the function f; (z) 
1 - 
Is=z 5 J G(z, C: d)f.(C) dl, 


and estimate the difference f (z) — f; (z). 
According to CaucHy’s formula 


fO=s5 1 GLAS OM. 


Therefore 


SA) h@=sq J CE.LIS Od? 


1 
toa G(z,C,d)[F()-f(O)] de. (2.6) 


In virtue of (2.1) and the second property of CARLEMAN functions the first 
integral in the right-hand side of (2.6) satisfies the inequality 


IJ Gf (Qdl[|<6-M. (2.7) 


In virtue of (2.2) the second integral in the right-hand side of (2.6) satisfies 
the inequality 
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I LF O-LO] Al sn, d)-e-7’, (2.8) 
where y’ is the length of the curve J”, 
1t(z,6)=max |G]. 


CeD 
Substituting (2.7), (2.8) in (2.6) we obtain 


I ; 
If (z)-—f5(2)| <5 [5M + 1 (z,d)-e-y']. (2.9) 
Let the CARLEMAN function be known. We denote the root of the equation 


6 ey 


ie) M ae 
by 6 (z, €), and put 
d=6(z,&) 
in (2.9). Then the inequality (2.9) becomes 
1 
If (2) — fs) S—M 4 (z, é) (2.11) 


It is evident that as e — 0, 6 (z, €) and therefore the right-hand side of (2.11) 
tends to zero. Hence, if the CARLEMAN function G (z, ¢, 6) is known, the func- 
tion f;, (z) may be considered to be an approximate solution of the problem. 
The difference between the ‘“‘true’’ solution f(z) and the approximate one 
fs, (Z) is thus estimated by the above inequality (2.11). 

It is evident that the smaller the function y (z, 6) is, the higher will be the 
precision of the approximate solution obtained by the above method. We shall 
now display a family of CARLEMAN functions with minimal y (z, 6) in the case 
where D is a simply connected domain. 

Let 0, (z, &) be the function which is harmonic with respect to the variable ¢, 
regular in the domain D and on J’ takes the boundary values 


0, (z,¢)=1n |z—€I, Cer. 


We denote by 0, (z, ¢) the harmonic function conjugate with respect to 
0, (z, 6) and by 6 (z, ¢) we denote the analytic function 


0(z, 0) =, (z, 0) +102 (z, 0 — 9, (z, z) — 10, (z, z). 


Let w (z) the harmonic function conjugate with respect to the harmonic 
measure @ (z), and W (z) the analytic function 


W (z)=w(z)+i@(z). 


2 Springer Tracts, Vol. 11, Lavrentiev 
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We consider now the function 


1 i" (z,z)+In6 


G60) tP 1—a(z) W400, -a df 


where 


a Gz)+nd-0), [AGN ge—40,¢,2)} 


g2,8)= i—o(z) 1—w(z) 
It can be easily seen that this function is a CARLEMAN function, 
0, (z,z)+1n6-@(z) 
= cp eee, 2.12 
u(z, 0) exp | 1—a@(z) ( ) 


The constructed CARLEMAN function is best possible from the standpoint 
of estimating the precision of the approximate solution by (2.11). However, 
it turns out that for certain domains the construction of the function is some- 
what cumbersome. Therefore, in order to solve the problem it is preferable, 
in most cases, to use another CARLEMAN function, namely ,the one pointed 
out by CARLEMAN himself (see. ref. [3], [29]). This function does not give high 
accuracy, but its actual construction is rather simple. 


§ 2. The Cauchy Problem for the Laplace Equation 


The Caucuy problem for the LAPLACE equation in the plane is similar to 
that of analytic continuation considered in § 1. An analytic function is te be 
determined from its values on a curve on which Caucuy data are given. 

In fact, let the value of the harmonic function u(x, y) and its normal 
derivative 0/On u (x, y) be known on some curve J’. We denote by f (z) (z = x 
+ iy) the function f(z) = u + iv where v is the function conjugate to wu. It is 
well known that on the curve J" 


v(z)= f Fu (2)ds+C 


where 2, is one of the endpoints of J’. 

Hence, if u (z), 0/On u (z) are known on J’, one may consider that the 
values of the analytic function f (z) on I’ are known. 

This section deals with the extension of the CARLEMAN function method for 
the case of harmonic functions in a thre-dimensional space. 

Let D be some domain of three-dimensional space bounded by a sufficiently 
smooth and closed surface 2. Let 2” be a part of X, and 


Ags pe Pe 
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We call a CARLEMAN function for the domain D and the surface Y any 
function of the three-dimensional vectors x, & and the scalar parameter 6 
satisfying the following conditions 


1) G(x, €,6)=—-—. + G, (x, €, 6) 


aS ,¢) 
where r (x, §) is the length of the vector x — é, and G;, (x, é, 6) is a harmonic 
function of €, regular and bounded together with its gradient in the domain D. 


2) The function G (x, &, 0) satisfies the inequality 


1 
2 feof 


where 0/On; is the normal derivative with respect to the surface 2; and do; is 
the element of area on & with respect to the variable &. 

Let the CARLEMAN function G (x, &, 6) for the domain D and the surface ¥ 
be known. By making use of the function G (x, £, 6) we can prove a “‘stability”’ 
theorem and carry out the method of concrete solutions of the Caucny problem 
for the LAPLACE equation in three-dimensional spacc. We dcnote by 4 (x, 0) 
the function 


Theorem. Let the harmonic function u(x) be regular inside D and contin- 
uously differentiable over the closure of D, and let it satisfy the inequalities 


~¢ 


\ does (2.13) 


ju(x, 0) = max fici+, + an G 


get’ 


lef w Ox) <e xe’ 
“ (2.14) 
wcol+}2weo <M xed” 
Then the inequality 
Ju (x)| <2 Mt(x,&) (2.15) 
also holds; where t (x, €) is a root of the equation 
u(x, t)-e=Mt (2.16) 


It is evident that t (x, €) > Oas e- 0. 

It is clear that the value of the harmonic function uw (x) at a point x inside 
D is expressible in terms of the uv, 0/Onu values on the boundary D accor- 
ding to GREEN’s formula 


1 0 0 As 0 ) 
u(x)=z— J \ogu-Zo- uh dogg J \¢- - u—z—G- updo, 
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from which, by virtue of (2.13) and (2.14) it follows that 
ju(x)|<u(x,d):e+M-6. (2.17) 


Assuming that 6 in (2.17) is equal to a root of equation (2.16) we obtain 
inequality (2.15). 

The construction of the concrete solution with the aid of the CARLEMAN 
function is completely similar to that given in Section 1. Let the Caucny data 
for the function w (x) on the surface X” be known with accuracy e, i.e., the 
functions @, (x), y, (x) are known on the surface 2” so that 


1p, (x) -—u(x)| Se, 
xed’ (2.18) 


<6, 


valu (a) 


and let the function (x) on the surface 2” be known to be bounded together 
with its normal derivative so that 


|e (x)| + 


u(x) <M, xeL (2.19) 


We consider the function 
us()=z- J | G-veO-26-9.(6 |do 
6 =] Ts! We An QP. é9 
and estimate the difference u(x) — u; (x). Substituting the corresponding 


integral expressions for wu, us and utilizing inequalities (2.13), (2.18), and (2.19), 
we have 


1 re) 6) 
\u (x) — ug (x)| = AS E at an G-u| do. 


(2.20) 


S 
<2M6b-+u(x,d)e-7— 


where S’ is the area of the surface 2”. 
Let 6 be a root of equation (2.16). Then inequality (2.20) takes the form 


4x 


Hence, if the CARLEMAN function G is known, the function u; (x) will be an 
approximate solution with accuracy estimated by the inequality (2.21). 


\u (x) —us(x)| < Mt (x, €) (242), (2.21) 
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We now tum to the problem of constructing the CARLEMAN function. The 
existence of the CARLEMAN function for an arbitrary simply connected domain 
D and for any part of the boundary 2” bounded by a smooth curve follows 
immediately from results of MERGELYAN, [47]. In that work it was established 
that for any arbitrary pair of continuous functions @ (x), wp (x) defined on a 
smooth piece of the surface X” and for any 6 > O there exists a harmonic 
polynomial H,, (x) of order m satisfying the inequalitics 


|H,(x)—@ (x)| <6 


Hal) v9) (2.22) 


<6. 


In this paper an algorithm is given for the construction of the poly- 
nomial. 

We consider the function [r (x, £)}-!, and let the point x lie inside thc 
domain D. We fix 6 and construct a harmonic polynomial P,,, (x, €) of order 
with respect to the variables ¢, satisfying the inequalities 


ame 6 C)— Pa (X; é)| <4 
bee" 


a1 . 0 
an! (x, 5) ~ a> P(X, o) Sd. 


We denote by G (x, &, 6) the function 


G(x, &,d)=r7*(x, €)— Pm (x, €). (2.23) 


It is evident that the function G defined by equality (2.23) is the CARLEMAN 
function for domain D and surface 2”. 

The above method for the construction of CARLEMAN functions is distin- 
guished by its great generality. However, the algorithm for such a construction 
is very complicated. We now give a simpler method for the construction of a 
CARLEMAN function for certain values of x. 

We denote by 2 (x, 7) the sphere with center at the point 7 and intersecting 
the point x. Let the point x€ D satisfy the condition: There exists a point 77 
such that the surface 2” lies inside the sphere 2 (x, 7). By expanding the 
function r~! (x, €) in TAYLOR series in the variables ¢, — 7,, we obtain 


Hoss hontsn (2.24) 


where P, is a homogeneous harmonic polynomial of order k in the variables 
&; — n; (J = 1, 2, 3). We denote by G,, the function 


G, (x, 4, €)=r7* (x, €)- d, Pa(%sme—m). 
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It may be casily seen that the function G,, satisfies the inequality 


|G,| + grad G,| < an |e (2.25) 
x 


1 2° 
|x—nl"** (Ix—nl-lé—nl) 
In fact, on any ray € — 7 = 2 ty where 7 is a scalar and y is a unit vector, the 
sum 


oy P,.(x,4,€-)= >. P,.(x,,T, ¥) 
0 0 


is a partial sum of the TAYLOR series in the variable t for the analytic function 
1 


F(a,x,4.y)=—_—_—_—__——"——. 
a a ae ae Sr 


The inequality (2.25) then follows from well know estimates for the reminder 
term for elementary analytic functions. 
From (2.25) it follows that if the number n satisfies the inequality 


9 n u" 
es <6 (2.26) 
Ix—yl""" [x-n|-R(Z)] 
where S$” is the area of the surface 2” and 
R(2")=max |E—y], Cex. (2.27) 


then the function G, is the CARLEMAN function for the domain D and the 
surface 2”. 


§ 3. Determination of an Analytic Function from its Values on a Set Inside 
the Domain of Regularity 


Let f(z) be an analytic function, regular in the unit disc D, 
If(z)|<1, zeED (2.28) 
and let A be some set inside the disc 
lzZ|<R<I (Dp). 


The problem is to determine f (z) in the entire disc Dp from its values on A. 
It is evident that for the solution of our problem it is sufficient to determine 
J (2) on the circumference 


|z|=R (FR) 


since according to the Caucuy theorem 
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Fac. (2.29) 


We denote by £ the operator which transforms the function given on I", 
to the function given on the set A by the formula (2.29), and by L, the operator 
which transforms, by means of the formula 


I Ja5— Wee Fat Wade, 
My 
the function given on J’; to the function given on Jp. It is evident that L, L, 
are linear, completely continuous operators. 
Solving this problem is equivalent to solving 


Lo=f (2.30) 


where f is the known function defined on A, and¢g is the function to be found 
on Jp. and 


p=Lyyp (2.31) 


where yp is a function defined on J’; satisfying the inequality 


Ipl<1. (2.32) 


Hence, the solution of the problem is reduced to the solution of an integral 
equation of the first kind (2.30), i.e., to the determination of the function o 
from the known function f, provided it is known that the function ¢ satisfies 
the conditions (2.31), (2.32). 

Now, we prove a theorem which characterizes the stability of the solution 
of the problem. 


We introduce first some definition sand notation. Let A be a set of points on 
some curve. Let us denote by H, (A) a system of 7 intervals containing the 
set A 

H,(A)>A 


and by yz (H,) the linear measure of H,. We call the infimum of  (H,) for all 
possible H,, containing A the n-measure of the set A and denote the 7-measure 
of A by 

ut, (A) =inf (A, (A)). 


It should be noted that this m-measure is unusual since, in general, it does 
not possess the property of additivity. However, it is evident that n-measure 
possesses the following property which replaces in some sense the additivity 


lay (Ay) + Lh, (Ax) > Ln, +05 (A, + A>) : 
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Now, let A be a set of points of the complex plane lying in the disc Dp. 
We denote by 4, the radial projection of the set A onto the circumference 


|jzj=r. 
Let W(z) be a conformal mapping carrying the unit disc onto itself, and 
carrying the set A onto the set A, c D,, where D,, is the annulus 
0<r<|zi<e<1 (D,o) 


By ul? (A) we denote the supremum of 4, (A,,) taken over all possible 
mappings W (z) 


pn (A)=sup Ly (Ayr) « 
Theorem. Let A be a set lying inside the disc 
lzi<R<1 (Dr) 


and let the analytic function f (z), regular in D, satisfy the inequalities 


If(l<l (zeD) 
If (ise (ze A) 


Then the inequality 


Inc(1—|z 
\f (Zl <exp {-c,med=l2)| (2.33) 
In u?? (A) 
is valid, where C, is some constant dependent of the radii 9, r, R and p satisfies 
the inequalities 
Hy (A) PP | Mp1 (A) 
< 
| C, <ex< C, (2.34) 


with C, also dependent on a, r, R. 
First we prove two lemmas. 


Lemma 1. For any points z, a lying in the annulus D,, the following ine- 
quality is valid 


a-Zz a,—Z, 


l—az| |1l—a,z, 


where a,, z, are the radial projections of the points a, z on the circle I’, 
FA lees ar 


Without loss of generality we may consider the number a to be real and positive. 
Let first |z| > a. We consider the function 


a—tz |? 
1—atz 


I,(t)= 
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Differentiating the function /, (t) we obtain 


Kas 5 —a’*)[(x* + y?)(1+a’? —ax)—ax] 


|1—a-z|* (Z=x+1y) 
Let 
x<a 
then 
(x*+y?)(L+a?—ax)—ax>x? 4+ y?—a?>0. 
Now let 
x>a 
then 
(x7+ 37)(L+a*—ax)—ax>(x—a)(1—ax)x>0. 
Hence 


IT, (1)>0 


from which it follows that for |z| > a, 


a-Zz a Za] 
a ve pond 118 2.35 
1l—az| |1—az), ( ) 
Now, let 
[z| =|a]. 
We consider the function 
2 
a—Z 
is 
ET eae —Taz 
Differentiating the function J, (t) we obtain 
Leo +y7]1-a%) 
1) aaa: e+ 
(i- 
from which it follows that for |z| = |a| 
a-Z a,—Z 
———| >|]. 2.36 
1—Gz| |1—4a,z, ( ) 


Finally, for |z]|< |a|, we consider the function 


2 
T—Z 


I,(t)= 


1—tz 


Differentiating the function J5 (t) we have 


,_ la? (x2+y7)—2axt 1] (a—a) (0? +y?—2axta*)(a-lel?—x) 
ie |1—az|* 
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Let 
a(x?+y*)—x<0 
then, clearly 


I,(a)>0. 
Now, Iet 
a(x*+y7)—x>0 
then 
ee: 
a—x>a(x°+y*)—-x 
and, hence 


Dino 2 Bae 2 
—x1A1—-— cs 
Fe at i A cs ER 

\l1—az*| 
from which it follows that for |a] > |z| 


a—z a\,;—2 


(2.37) 


l1-—az 1-4a),,°z 


The inequality to be proved follows from inequalities (2.35), (2.36), and 
(2.37). 


Lemma 2. Let a,, az, Z,, Z, be some complex numbers, and 
|a,|=|a,|=|z,;]=|z.|=r<1. 
Then, if the inequality 


arg—-<arg—*<x 


ay a 
Z1 Z2 


is true, then the inequality 


will also be true. In proving the Lemma, one may suppose that 
a,;=a,=a 
where a is real and positive. 
Let 
z,=ae'*!; z,=ael®. 
We consider the function 


Z 2 


1 


Differentiating the function J (@) we obtain 


I(g)= z=ae’’ 
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; a’ sing (1—a’?)? 
I ae ee 


from which we obtain the inequality to be proved. 


We now prove the theorem. Let w = w (z) be the conformal mapping of D 
onto itself for which 


Ly (A) = Hy (Ayr) 
Under this mapping, the set A is carried into the sct 
Ay ¢ De 


Let us consider the function 
o(z)=In|f (z)|. (2.38) 


It is known that @ (z) is a harmonic function, regular at those points of D 
where the function f (z) is not zero. We consider some point a belonging to A,,. 
According to the condition of the theorem 


|f (a)| =e,<€ (2.39) 


From (2.38), (2.39) and well known properties of harmonic functions it follows 
that 


If Ql=ea 


on some analytic curve y, intersecting the point a (if f (a) = 0 then the curve 
consists of the single point a). 
Let the curve y, be an open arc in the annulus 


te 


2 1 » 
—r<|z]<— 2 D 
z 7S Is3 (+ 0) 


Then the curve y,, either intersect one of the two circles 
\z| =a (F,) 


1 
l=z+20) (F) 
or bounds a domain containing D2jar- 


Suppose 7, intersects the circle I’, in some point a,. The distance between 
the points a, a, satisfies the inequality 


la-aler. (2.40) 
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It is not difficult to show that by virtue of (2.4) the harmonic measure of the 
curve y, in the domain D with respect to the point z satisfies the inequality 


r. S=|z| 
aa ally 2.41 
Aas OES Pa aa 
From (2.41) we obtain 
PO iseotne 2 (2.42) 
a ee care 


In case the curve y, intersects the ring 
j 
|z|=z(1+2e) 


we similarly obtain the inequality 


1-—@ 1-|2| 
<exp \Ine—_— ———>. 2.43 
If ise P| Sra (2.43) 
Finally, if the curve y, bounds a domain containing D2/3, then from the 
maximum modulus principle, in the whole domain D2/3,7 the inequality 


If@l<e zeD2, 


is valid, from which the assertion of the theorem immediately follows. 

Thus, if just one of the curves y, is not closed in the annulus Dio the 
statement of the theorem follows. Now suppose all the curves y, are closed in 
Dy» and consider any curve y,. As y, is closed in D, it follows that inside the 
domain bounded by the curve y, there is at least one zero of the function f (z). 
We denote by y the sum of sets y,, and by y, the projection of y on the circle 


an 
le, 


|z|= 


In the ring 


i = 


te 


suppose the function f(z) has p zeros at the point a, (k = 1, ..., p) (roots are 
counted according to their multiplicity), and let 


F@=f@) IL 


—A,Z 
a,—Z : 
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We prove now that in D,, there exists a point z* such that 


| p 
F (2* ——__ 
[F ( 1s (2.44) 


Let z be some point of Digs and z, be the radial projection z onto F 
By virtue of Lemma 1 we have the inequality 


(2.45) 


Let us consider now the set y,. Since the function f (z) has p zeros in Dis 
it follows that the set y, consists of no more than p intervals. Then by virtue of 


Wr? Av, 
it follows that the measure of y, satisfies the inequality 
(y= My (A). 
We denote byy; the part of y, which belongs to the semi-circle 0 < argz< 2, 


and by yy the difference 


t 


Ye =Vr—Y; 
It is evident that one of the two inequalities holds 
ioctl, ise. 
HO)E> HO) 5 


um” 


1 
i (yy) = ox a (y,) 5 
For the sake of definiteness assume 
; 1 
LO)DZ> HO) - 
By virtue of Lemma 2, in proving (2.44) one may consider y, to consist of 


one interval. From CHEBYSHEFF’s results on polynomials which deviate least 
from zero it follows that on the interval y} there is a point z¥ such that 


[leer — 241 [Cle Gy ]P= [euce (ayy, (2.46) 


where C, is some constant which depends on r. From (2.46) and (2.45) we see 
that on the set y we can find a point z* such that 
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£ 1 —a,z* 


i ; (2.47) 


| 1 
LC, (A) 
and the inequality to be proved follows from (2.44). 


Now we consider the function 


v(z)=In |F (z)| 


1 a,—z* 


The function F (z) has no zeros in the annulus D,» and therefore the function 
v(z) is a regular harmonic function in D,,. Then by virtue of (2.44) and the 
conditions of the theorem 


v(z*)< —pInyp(A)+Ine 


v(z)<O = (ze D,,,) Ce) 


Wc denote by Fro the boundary of the annulus D ro: It follows from (2.48) 
that 
J v@ds<C,, [Ine— pln p? (A)] (2.49) 
Te 
where C,, is a known constant dependent on the radii r, @ and it follows from 


(2.49) that the function v (z), in the annulus D,,, satisfies the inequality 


v(z)<C,, [Ine— pin pi (A)] (2.50) 
From the inequality (2.50) we obtain 
E Cre 
lf Dis) ——,} (2.51) 
Lup’ (AD)? 


It is not difficult to obtain from the inequality (2.51) the following estimate for 
the modulus | f (z)| in the whole disc D 


FI } ee 
7) | ES eee (2.52) 
Lup (A)]? 
Then, since the function f (z) has p zeros in the annulus D,, it follows that 
If (Dlsv’, (2.53) 


where 
_ll+8. . 2+ 


) = —— 


ieee 83 


~ 


L 


It is easy to see that 
Inv<C,-(1—|z|) (2.54) 


where C, is a constant dependent on a. 
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Let p satisfy (2.34). Then by virtue of (2.52), (2.53), and (2.54) the inequality 


Ine 
lf @sexp}—C (1 -ep—e_ 
In py (A) 
is valid, and this concludes the proof. 


We cite some examples of obtaining pre (A) and the corresponding esti- 
mates. 


1. Let the set A,,, have positive JORDAN measure # (A,,,). Then, it is evident 
that 


ey (A) 2 Ht (A mh, ; 
it, (A) > Hu (A,,,) : 


n> 0 


If, in particular, A,,,, consists of curves, then 


Uy (A)> U(Av,) » n<m, 
iy (A) = E(Avr), n>=m. 
If m1 = 1, then the estimate of the theorem, accurate up to the constant C,,. 


coincides with the estimates for an analytic function using harmonic measure. 
2. Let the set A consist of the points a, (k = 1,..., ) 


la,|=a. 


arga,=ko, (s==). 


One may easily see that in this case 
L,(A,)=a(n—p)o, psn, 
H,(A,)=0, p>n. 
The relations (2.34) in this case are as follows 


[(n—p)d-a]’<e°*<[(n—p+1)d-a]? (2.55) 


It is not difficult to show that the number p satisfying (2.55) satisfies the ine- 
quality 
1, nd 
pec ne 
1+6'Ine 


from which we obtain 


nd‘ Ine’ 
< C,————(1--|z ; (2.56) 
lf (z)| P| ia snes \z|) 
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3. Let the set A consist of the points a, (k = 1,..., 00) 


lal=a, argay=—. 
2 
One can see in this case that 
a ; 
Hl A)=—, (2.57) 


and the relations (2.34) are as follows 


from which 


p*>Cz,/[inel 


and 


If (2) <exp{—C, J/ In] el-(1—|z1)} (2.58) 


§ 4. Analytic Continuation of a Function of Two Real Variables 


Let f (x, y) be a function of the two real variables, analytic and regular in 
the disc 


x*+y?<1, (D) (2.59) 


and let A be some set inside D. Let the values of the function f (x, y) be known 
on the set A and suppose it is required to determine f (x, y) inside D. We 
prove a theorem characterizing the stability of the solution for a certain class 
of sets A. 

First some definitions are introduced. Let us divide the x, y-plane into a 
system of parallel strips H, of width 4, (k = +1,..., + 0). Let T denote this 
division while A, denotes the part of the set A that belongs to H, 


A,=ANH, 


By A, we denote the projection of the set A, onto the center line of the strip 


H,. Let t> 0, C;(j = 1,...) be some constants and let , be an integer 
satisfying the inequalities 


aA ne t1 7 Nk 
E =s(4d)] <t+ Ch, < ES (2.60) 


where 44, (A,) is the n-measure of the set A,. 
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We denote by H,, (A, T) a set of strips A’, for which the following inequality 


Is true 


In(t+C,h,) 
sl Si Aid LA 2.61 
CyInp,, (Ade een 


and by #,,,, (A, T) the m-measure of the projection of the set H ar (A, T) onto 
some straight line perpendicular to the lines of the division T. We call the supre- 
mum of #,,,, (A, T), for all possible divisions 7, the (m, 2, t)-measure of the 
set A and we denote it by 


Linde ssp {imiae (A, T)} . 


Theorem. Let A be a set of points of the x, y-plane lying in the disc 
x?+y7<R?<1 (Dz) (2.62) 


and let f (z,, z,) be an analytic function of two complex variables z,, z, regular 
in the hypercylinder 


(Rez,)*+(Rez,)*<1, 
lImz,|<x, |Imz,]<x. oe 


If the function f (z,, z,) satisfies the inequalities 


Lf (21,221 <1, (Z1,Z,)EQ, 


(2.63) 
If yy] <e (G,yed, 
then the inequality 
Lf (x, y)|-<ex =i (x, y)eD (2.64) 
vy => p In Anre(A) ’ oy R ar 
is valid, where the numbers m, A satisfy the relations 
m+i1 m 
Un+1 A, (A) <2? < Linas (A) (2.65) 
C, 2 
and C; are constants depending on R, x. 
We denote by 7, the division of the plane x, y for which 
Imre (A) — Ponae (A, T)| <a. (2.66) 


One can obviously consider the corresponding strips H;, to be parallel to 
the x-axis. Let the center line of the strip H, be the straight line 


Y=YVx 
It is not difficult to show that from (2.63) it follows that 


3 Springer Tracts, Vol. 11, Lavrentiev 
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lgrad f (x, yI<C3, (x, y)eDpR (2.67) 


where C3- depends on R and x. From (2.67) and (2.63) we find that at the points 
of the x, y-plane belonging to the set A, the function f (x, y) satisfies the ine- 
quality 


1 
If yisets Cahy- (2.68) 


Let us consider the function 
Py (Z)=Sf (2, Ve)» (x, y)EH,. 


From the condition of the theorem and from (2.68) it follows that the function 
@, (z) is an analytic function, regular in the domain 


[Im z| <x; (2%) 
and satisfies the inequalities 
lp, (z)| <1, zEQh; 2.69) 
lo(z)l<et+Czh,, (Rez, y,)€A,. e 


Applying the results of the preceding section to the function @ , (z) we find that 
the following inequality is valid 


[Re z| <——— 


In(e+C3h,)-C, 
: < Se 
Ip, (z)| Scexp | In Lng (A,) ) 
x“ 
Im z|<— 
[Im z|< + 
where the constants C3 and C, depend on R, x hence 


In(e+C3h,)°C,4 
In Ln, (A,) 

Now we consider the function 

W.(z)=f (x, z). 


From the conditions of the theorem, the definition of the set H,,(A, T) and 
(2.68) it follows that the function w, (z) is analytic, regular in the domain 


If Ce lsexe | \ (x, y)eH,. Dp (2.70) 


[Re z|<,/1—x? 
(Q5) 
|Imz|<x 
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and satisfies the following inequalities 
Iv. (ZIS1, zEQ,, 
(2.71) 
lW.Qi<e*, (x, VY)EH,, (A, T,) 9 Dr. 
Applying, as previously, the results of the preceding section we obtain, by 
virtue of (2.71) 
(xv, Rez)eEDp 


ACz 
A(z) <exp 4 
x 
Im z|<— 
\Imz|< 5 
from which it follows that 
Pim miewot = (x. p)ED (2.72) 
oa : ID {nae (A, 1) , : . : 


Passing to the limit at « > O in the inequality (2.72) we obtain the inequality 
(2.64) to be proved. 

Now we give two examples of the applicability of the theorem. 

1) Let A be a closed set with positive plane measure sc (A). We show that 
in this case when 


a [Int|-C 
~ [In 2 (A)| 
the inequality 
1 
Hinde (Aa cH (A) (2.73) 
is valid, and consequently 
Ing 
3 < ey gre IO 2.74 
I ( vIsexp| oa (2.74) 


where C Is some constant. 
Let N be some number. We perform a division T of the plane into strips 


Ktsysx 
and denote by H,,; the square 
k= 
N ~"—N 
(H;,;) 
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We denote by p, the number of squares H,, of the strip H, containing 
points of A, and by p the total number of squares containing points of A 


Nn 
p= » Py - 
— mn 
One can easily see that the inequalities 
pp A): N* 
(2.75) 
—2m 


| Px 
——E 
Bin (A,) = N 


are valid. In fact, the first of the inequalities (2.75) follows from the definition 
of plane measure of the closed sets. In order to obtain the second of the ine- 
qualities (2.75) we note that by virtue of the definitions of p, and A,, on the 
straight line 


there are p, non-overlapping intervals of length 1/N containing points of A,. 
Thus the sum of the lengths of any intervals covering A, cannot be less than 
the number 
Pp-2m 
N 


from which the second inequality follows. 
Let us consider the strips H, for which the inequality 


1 
Pr2q N-H(A) (2.76) 


is valid. 
Let (2.76) hold for 


k=k, (j=1,...,q) 


It is not difficult to show that the number q satisfies the inequality 


a2 N- H(A). (2.77) 


In fact, assuming the opposite we see that 
1 2 
p<2Nq+(2 N—-q)' 7 N-n(A)< N*+ (A) 


which contradicts (2.75). 
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We denote by H the set of center lines of the strips H,; (j = !,...,q), and 
by H the projection of H into the straight line x = 0. 

In analogy with the second of the inequalities (2.75) we find that for any x, 
the n-measure of H satisfies the inequality 


q—2n 


A) > 
i, (H) > N 


from which, by virtue of (2.75) and the definition of s¢,,,, (A, 7), we obtain 
n< Eat N-p(A) 
‘6 


Haac(As TE HCA) (2.78) 


ne Int-C 
“In w(A) 


Passing to the limit as N — oo in (2.78) we obtain the inequality (2.64). 


2) Let the set A consist of points p;; with coordinates 


1 . 
ViGg=AT> i?) Peete’ 6) 

24 
1 = 
ae j=i,itl,...,00 


We show that in this case for A < Cy Vint the following inequality is valid 


1 
Minar (A) a am ET acer (2.79) 


and consequently, 


Lf (x, yy <exp {—C [In e|/*} (2.80) 
Let us perform the division T of the plane x, y into the strips H, 


k-1 k 
or Phy 
where N is some number. 
It is easy to see that for k = 2? (p = J,..., N — 1) the m-measure of the 


set 4, is equal to 


Ly (A= (2.81) 


gntp’ 


Therefore the relation (2.60) in this case are as follows 
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n+ Nk 
=m <n <( C1 ) (2.82) 
gm tert aie! P 


where 


From (2.82) it follows that for 


p<J|int,|, 2<C3,/|In7| 
the straight line 


k—0.5 
aN 


~ 


VaIn= (k =2") 


belongs to the set A,, (A, T). 
By virtue of the above statement, one can show that the measure of 
the projection of H,, (A, T) onto the straight line x = 0 satisfies the inequality 


1 1 
H,.(A, T)|>——- 
Lay [ acl 1 qm QV [in tif 
from which follows 
1 1 
A, T > - 2.83 
fonar ( ) am QV jin fil ( ) 


Passing to the limit as N — 00 in (2.83) we obtain the inequality (2.79). 


§ 5. Analytic Continuation of Harmonic Functions from a Circle 


Let u (z) be a harmonic function regular and bounded in the disc 
|z|<1 (D) 
and let its value be known on the circle 


Izj=r<1 (r,), 


u(re?)=f (9) 


where f (@~) is a known function. It is required to determine the function u (z) 
inside the disc |z| < r from the values of the function on I’,. Let us prove a 
theorem characterizing the stability of the solution. 
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Theorem 1. Let the considered function « (z) satisfy the inequalities 


|. -% , 
= J wre) dose, 


(2.84) 
nee re 
Br 2 el? <1. 
ju (e"")do<l] 
Then the following inequality is valid 
| eee oe Ine:Ing 
— | u“ (ee) dg <exp ;——— 
me (oe) de Py ia (2.85) 
It is well known that 
u (gel?) = ys o* (a, coskg+b, sin kg) +2. 
1 
Let us write 
=z 5 a 
Jagt be =Cy. [B= 
From (2.84) we have 
x 
y aud Ore 
Z (2.86) 
eal 
0 


Let us estimate the function 


w= Ye Ci=— J u*(0e) do 
0 -_ 


! 
t 
using (2.86). 

By virtue of (2.86) the sum to be estimated is bounded. Let the sum reach 

its conditional maximum over the variables C, for C, = C, (k = 0,..., 00). 


It is evident that 


C,=0, K+D,q, 


where p, 4, Pp < q, are some numbers, since otherwise it is quite clear that for 
C,, one may choose variations 6, such that 


» (C, +6, 07% > y C207 
0 0 


2 (Cito) < y Cys 


y, (CE, +5) 1 < YY Cer. 
0 0 
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Therefore C pe C, satisfy one of the three relations 


Ch+Cr=1, ko 


Cyr Cre: 


(2.87) 
C,=0; (II) 
C,=0. (II) 
Let the first of the relations (2.87) hold. Then 
ae e—r74 ; 
Pp peep? 
2.88 
ze 7 p2P_g ( ) 
i = peP = re 


We note that from (2.88) follows 
rt<cecr??, (2.89) 


By virtue of (2.87), (2.88) we obtain 


oo 2p 2q 2p,24q_ 24 ,2p 
m2.2K 0 =e r“@-—r “@ 
DCEO 8 aegypti, - (2.90) 
It is easy to show that for p < q the functions 
2p 2q 
Q-—?Pp 
F(R. Dao aa? 


2p,2 2 
r Po q__ > 19?P 


F2(p,q)= 


2 
per p24 


increase with increasing p as well as with decreasing g from which, by virtue 
of (2.89) and (2.90), it follows that 


sce nee : Ine-In 
Y Cres lim {eF (pg) +F 2 (p,a)} =exp jnetnel (2.91) 
pqrat. 
Inr 


Now let the second of the relations (2.87) hold. Then 
oO 
y C2 9?*— C2924 
. q 
0 
and by virtuc of (2.86) 
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from which we obtain 
Y. Cro**<exp {iseinel (2.92) 
Q 


In case the third of the relations (2.87) holds, incquality (2.92) can be ob- 
tained in a similar manner. 

The statement of the theorem follows from (2.91), (2.92) and from the 
definition of the numbers C,. 

With respect to the function u (z) considered in theorem 1, let it be known 
that at the boundary of the domain the n-th derivative is bounded in L,, i.e., 


—  § ie u | dp<M 


and suppose it is required to determine the function wu (z) on the circle 0 = | 
from the values of the function «(z) on the circle @ = r. We shall prove a 
theorem characterizing the stability of the solution. 


Theorem 2. If the function w (z) satisfies the relation 


fu? (re'*) dose; 


| dp<i; (2.93) 


J u(e’)dp=0, 


nt 


then the following inequality is valid 


Ls F 1 
ae i 7(e °) dos => qn’ (2.94) 
where 2 is a root of the equation 
2Alnr—2nIniA=lIne. (2.95) 


We note that for a sufficiently small ¢ the root of equation (2.95) satisfies the 
inequality 
1 Ine 
>— —. 2.96 
4 Inr ( ) 
The proof of theorem 2 is similar to that of theorem 1. From (2.93) it follows 
that 
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CO 
iy r7*C? <e; 
0 


: (2.97) 
y k*Cr<t, 
(0) 


n— 


where C,, as in the previous theorem, are expressed by the expansion coeffi- 
cients of the function u (ge'”) in a Fourier series in the polar angle @. 
In analogy with (2.87) we find that the numbers C, for which the sum 


5 Chat i u*(e'*) do 
fr) Tl —x 
attains its conditional maximum, satisfy the relations 
C,=0, K+D,q, 
while C,,, C, satisfy one of the three relations 
rPC2 +r°ICi =e ! C,=0 (I) 


pC +q"Cr=1 C,=0 (il) 


q 


from which the theorem follows in an analogous manner. 


§ 6. Analytic Continuation of Harmonic Functions with 
Cylindrical Symmetry 


Let x,y be vectors with components (%,,...,%,), (4. ---» Yert1) and let 
u(x, y) be a harmonic function regular and continuous in the domain 


Ixi<1,  Jylst @) (2.98) 


Let also the function u (x, y) as a function y depend on the variable r = |p| 
alone. It is required to determine the values of the function u (x, y) inside D 
from its values on the set P 


MVEP.. i. v=o: |x|<o<l. 


We shall prove a theorem characterizing the stability of the solution. 


Theorem. If the function u (x, y) to be considered satisfies the inequalities 


lu(x,yl<l, (x, y)eD; 
(2.99) 
\u(x, y)|<e, (x, ye P 
then the following inequality also holds 
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lu (x, y)|<Ce’, (x, yyED,, 


where D, is any domain lying strictly inside D and C, y are some constants 
ependent on D, and 0. 

By virtue of the conditions of the theorem, the function « (x, y) may be 
represented as a function u(x,r) where r = |y|. The function u(x, r) will 
satisfy equation 

O71 om Ou 


ape ty Bp TAstt=9- (2.100) 


From the known properties of harmonic functions it follows that « (x, r) 
is analytic with respect to the variables x, r evcrywhere inside D. 

We will estimate the modulus of the function A?« (x, 0). We consider the 
function f (z) where z is a vector with complex components z, = x, + if, 
(A = 1, ..., m) analytic with respect to the variables z, and coinciding on a real 
hyperplane with the function w (x, 0) 


f (x)=u(x,0). 


[t is evident that the function f (z) is regular and bounded in the domain R 
of the n-dimensional complex space 


zeR, if |Rez|<o, |Imz|<h, o<e,<1, 
where /: is a constant dependent on o@,. Define R, to be the domain 
zéER, if |Rez|<e, |Imz|<h,<h 


From (2.99) as well as from the results of § 1 it follows that the function 
f(z) in the domain R, satisfies the inequality 


If (z)|<Cye”; (2.101) 


where C;, y, are some constants dependent on 9, /, @,, 4, whereas from (2.101) 
and from known properties of analytic functions it follows that 


If (2) sk! ate” (2.102) 


where « is a constant. 
From (2.102) we obtain 


APu (x,0)<(2 p)! aie" (2.103) 
We now expand the function u (x, r) in a MACLAURIN series in the variable r 


ee ce (x,0)-1°? (2.104) 
uGnr)= ) Oarae. x,0)+7 os 


(the coefficients of the odd powers of r are zcro by virtue of the symmetry). 
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Substituting (2.104) in (2.100) we find that the functions 027/dr?? u (x, 0) 
satisfy the following relations 


07? 2p-1 or") 
Sap Saale aa (p=1,..., 0) 


from which 


4 


a 5 
a5 0=(—1) 
The validity of the statement of the theorem for any neighbourhood P 
follows from (2.104), (2.105) and (2.103). The validity of the assertion of the 
theorem for the arbitrary domain D, follows from the validity of this assertion 
for the neighbourhood P and from general estimates for analytic continuations. 


(2 p—1)!!(m—-1)!! 


Opti ee OD) 


Chapter Iii 
Inverse Problems for Differential Equations 


We consider two types of inverse problems for linear differential equations. 

The first kind is the so-called inverse problem of potential theory where the 
equation is non-homogeneous, the coefficients of the differential operator are 
given and it is required to determine the right-hand side. In this case we shall 
restrict Our consideration to the LAPLACE operator, although a number of 
results are carried over to more general elliptic equations. 

The second kind is that in which the equation is homogeneous and it is 
required to determine some of its coefficients. A well-known problem of this 
kind is the inverse problem for the STURM-LIOUVILLE equation. 


§ 1. The Inverse Problem for a Newtonian Potential 


By an inverse problem for a NEWTONIAN potential we mean the following. 
On a part of the surface of the unit sphere D, there is known a harmonic func- 
tion u (x) which is an external NEWTONIAN potential of a body with unit density 
which is star-shaped with respect to some point 7 and lies inside the sphere 

Ixj<R<1 (Dp) 

It is required to determine the dimensions and shape of the body from the 
function u (x). This problem has an important application in geophysics and 
a number of papers are devoted to its investigation. 

The uniqueness of the solution of the problem was proved by Novikov. 
The notion of correctness in a compact set was introduced by TYKHONOv in 
an application to a problem. In [16] estimates were obtained which characterize 
as properly posed the problem for a class of domains bounded by continuously 
differentiable surfaces with bounded derivatives. 

Let the polar coordinates of the vector y — 7 be denoted by 0, @, @ and let 
the surface of the body inducing the potential u (x) be defined by the equation 


e=f (@,9). 


It is well-known that the function u (x) is equal to 
n/2 mn f (9,9) o7 sin 0 


= —— dod dé 3.1 
wore LUT epee en 


where r (x, y) is the distance between the points x, y. 
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Equality (3.1) may be regarded as an integral equation for the unknown 
function f(g, 0) which defines the body to be found. If the integral operator 
in the right-hand side of (3.1) is denoted by A then (3.1) will take the form 


u(x)=Af (@,0). (3.1’) 


The operator A in the right-hand side of (3.1’) clearly transforms a function 
of the polar angles g, 0 to a function defined on some part of the unit sphere. 
Thus, the solution of the inverse problem of potential theory is equivalent to 
the solution of the nonlinear integral equation (3.1) or (3.1’). 

Now we proceed to the presentation of a stability theorem for the inverse 
problem of potential theory. For simplicity we restrict ourselves to the case in 
which the part of the surface D on which the potential « (x) is known, is defined 
by the inequalities. 

x, 2H, |H|<1 (S,) 


Theorem 1. Let the functions fi (@, 6), f2(@, 0) satisfy the following con- 
ditions: 

1) The functions f;(@, 0) (j = 1, 2) are continuously differentiable for all 
values of the polar angles g, 6 and their derivatives satisfy the inequalities 


lo a 
aq Jil? O)|<M; eae. 0)|<M (x*—67). (3.2) 

2) The following inequality holds 
[u, (x)—u,(x)|<e, xESy (3.3) 


where 
u;(x)=Af;(Q, Q) : 
Then the difference f, (y, 8) — f, (Q, 9) satisfies the inequality 


C 
If: (9, 8) — fa(@, 1 <— (3.4) 
where Cj is a constant depending on M, and n is a number defined by the 
relations 
1 2n+3 1/C3 1 2nt+1 
rats <C, }Ine <(+) (3.5) 


while C; are constants depending on R, M. 
We first prove a lemma. 


Lemma 1. Let w (y) be a summable function which does not exeed one in 
modulus and let the function f (@, 0) satisfy the conditions of the problem as 
well as the inequalities (3.2). Then if the function 
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nf2 xn f(,0) 


2. 
sind:@ 
u (x)= f f f Q (y) 
—n/2—-zx 0) r(x, y) 
satisfies the inequality 


do dg d0 


Ju(x)ise (veS,,) (3.6) 


then the following inequality also holds 
1 


ju (x)|+igrad u(x)]}<C, [Ine] <s (3.7) 


where the constant C; depends on R, M and x is any point lying outside or on 
the surface 


o=f (9,9). 
It may be shown from (3.6) that the function « (x) satisfies the inequality 
Ju(@ise%=e, ([x[21) (3.8) 


where Cg is a constant depending on the radii R and H. 
Let the polar coordinates of the vector x — 1 be denoted by 0’, g’, o’ and 
the distance between the points x and y by 


r(o', 9’, 0’, 0,0, I)=r(x, y). 


Let us consider the function 


n/2 « f(@,0) 2¢: 
, t QO sin 6-@(y) 
o(z,9',0)= f f 7 


——..—— dodgdd. 
—n/2—-n 0 r(Z,@ , 0’, 0, @, 9) ? 


It can be easily seen that w (z, 9’, 0’) is an analytic function of z which is 
regular in the domain 


Rez—C,|Imz|<f(@',0') (Qg4’) 


where C7 is a constant depending on the constant M in (3.2), and satisfies the 
following inequalities 


|o(z, o’, 0')|+|grad w(z, 9’, O)|< Cg, ZEQg 6 ; (3.9) 
d , , t , 
“grad (2, 9',0')}<ColIn[z—S(9', OV], aa 
Rez>f(qg’,0’); Imz=0; 
la(z, 9,0 )|<e, 
(3.11) 


Rez>1, Imz=0, 


where C; are constants depending on M. 
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From the inequalities (3.9), (3.11), some properties of analytic functions 
and the results of Chapter II, it follows that the function © (z, 9’, 0’) satisfies 
the inequality 

lo (z, 9’, 0) + |grad w(z, 9’, 0’)| <e°"* 
3.12 
Rez>1, Imz<—C,. ( ) 


By virtue of (3.9), (3.10), and (3.12) the functions  (z, g’, 6’) and grad w 
satisfy the conditions of Carleman’s theorem ((3] or [29]) and therefore the 
following inequality holds 


|o(z, 9’, 0’) + |grad w(z, 9’, O')|< Cy, exp {Ine |z—f (9’, 6’) “7}. (3.13) 
Using inequalities (3.10) and (3.13) we obtain 


[u(x)| +lgrad u (x)|=lo(e', 9’, 0+ 


+|grad w(9', 9’, OS Cy, exp {Ine[O—f (9', OI} + (3.14) 
+C6lo—f (9’,6)] {1+ llnLe—f (9, O)]1; 
f(g’, @)<e'<@<iltyp. 

Setting 


1 
0—f (9',@)=|Inel 2c: 
and using the obvious inequality 
1 
7 t or = t anya ©12 
lin [@—f (9, OVNI <C,3[@—-F (@', OT 
we get the inequality (3.7). 


Lemma 2. Let w (x) be a harmonic function regular inside the surface 
Ix—nl=f(@',0)+h, h>o (2) 


where f (@, 9) satisfies the conditions (3.2) and does not exceed one in modulus 
in this domain. Then on the surface 


Ix—nl=f (¢', 0) (Xo) 
the inequality 
070 dw] C, 
+¥l—|< 3.15 
d Ox OX, Lex, oy) 


holds, where the constant C, depends on M. 
It is evident that by virtue of (3.2) for any point x, € 2, there exists a sphere 
of radius C,h with center at the point x) lying completely inside the surface 2’, 
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(C, depends on M). By the maximum principle the values of w (x) at the 
boundary of this sphere do not exceed one in modulus. The assertion of the 
Lemma follows from the above and from known estimates for derivatives of 
harmonic functions at the center of a sphere of given radius in terms of its 
values at the boundary of the sphere. 


Lemma 3. Let w (x) be a harmonic function regular inside the surface 2, 
of Lemma 2, satisfying the following conditions. 


1) For some vector x) € 2, and for an arbitrary point x lying on the surface 
2, and inside the sphere 


|x —Xol <i, 
the function w (x) is 
wo(x)=1. 
2) For an arbitrary xé 2’, 
|w(x)} <1. 


Then the following inequality holds 
h \° 
1—w(Xo,)<C, (+) (3.16) 


where C; are some constants depending on the constant M; x9, — is the 
vector with the polar coordinates 


Co—h, Po, Oo 


where 0p, Qo, 99 are the polar coordinates of the vector x9. 

From the conditions imposed on the function f(@, 0) and also from the 
maximum principle it follows that the assertion of the Lemma is equivalent to 
the following assertion. 

Let W’ (x) be a harmonic function regular in the sphere with a cut-out cone 


o<A, Oza, 


(e,”, 9 are polar coordinates of x — and « is a constant depending on M). 
Suppose W’ (x) takes on the following boundary values 


W'(x)=-1 
on the surface of the sphere ge = A; 
W'(x)=1 


on the surface of the cone x/2 — 0 = «. Then at the point x,, where the polar 
coordinates of the vector x, — 7 are equal to 
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v1 


Q,=h, 9,=9, =z 


the function W’ (x) satisfies the inequality 
h \ 
1—-W'(x,)<Cy i > 


where C; are constants depending on «. 
The latter assertion follows from known estimates for harmonic functions. 


Lemma 4. Let w (f) be a nonincreasing function of the parameter ¢ on the 
interval [0, C] 
@a(0)=1, w(C)=0. (3.17) 


Then for an arbitrary t,0< ¢< 1,n > 1 thereexists anumberk4,1< k<n 
such that 


o(t!)—o()<—. (3.18) 


We consider the sum 


ps [o(t!*)— a(t] =o(f)-o(). 


From the conditions of the lemma all the terms of the sum are non-negative 
and the whole sum does not exceed unity, from which it follows that some 
member of the sum does not exceed I/n, i.e., for some k (3.18) holds. 

We now turn to the proof of the theorem. We introduce the following 
notations. 

Let 


f' (9, 0)= Neri LY; (9, 6)] ) 
f'(9,6)= min [F,(9,6)]. 


We denote by 2? the surface defined by the equation 


o=f'(9,4) 
and by 2? the surface defined by the equation 
e=f' (9,4); 


and by x ny the surface defined by the equation 
o=f'(9,0)+h. 
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We denote by zh the part of the surface ai) whose points satisfy the inc- 


quality 
fi (9, 0)> fr (Qg, 0); 

and by xin we denote the part of xin whose points satisfy 
ti (9, 0)<f, (9,0). 


The corresponding parts of the surfaces 2?, Y‘ will be denoted by 
yt yi- pit ie 
We denote by Di , the part of Zi whose distance from ais is not less than 4 
r(x, Sen XE Liha 


ee denote by py (ay a the part of ais which belongs neither to bares » nor to 
>the , and finally, by Zi) ay a, We denote the part of a ) Which does not belong 
to the surfaces 
Zoi az» Diy dar Zn ay (A > Ay). 


The parts of 2, 2? to which correspond rays from the point 7 intersecting 
the indicated surfaces, will be denoted by 


zy" a> Lp Li, A2 
zy a Zs Dae 
respectively. 
We denote by 2 the domain enclosed between the surfaces 2!, &* while 
by QF , 27 we denote the parts of 2 enclosed respectively between the surfaces 
Pe zis and £)-, Zi, and the radial straight lines bounding them. 
Finally we denote by u (2) the volume of the domain 2, by o (2) the area 
of the surface 2, by « (x) the function defined by the relations 


a(x)=1 xEQ5; 
a(x)=—1 xXED, 5 
a(x)=0 outside Q; 


and by u(x) we denote the difference 
u(x)=u, (x)—t2 (x). 


ao) 4 
r(x, as 


Let v (x) be a harmonic function regular inside the surface 2” and contin- 
uously differentiable on 2'. Then from potential theory 


Clearly 


u(x)= J 
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facpoepdxe ge f(SuoZo-u)do (3.19) 


Q Zi 
where 0/0n is the derivative along the normal to 2%, 
We set 
»(x)=div[W(x)(x—n)] 


where W (x) is a harmonic function regular inside xn) and assuming on on) 
the following boundary values 


W(x)=1 xe Xa; 
W(x)=-1 xEXGy4; 
W(x)=0 xEZiy4 


We estimate first the right-hand side of (3.19). By virtue of Lemma 1 and 
the conditions of the theorem, the following inequality holds on the surface 2” 


|u(x)|+|grad u(x)|<e,=C, fiivel ee: (3.20) 
and by virtue of Lemma 2 
[o(2)1 + Igrad 0(x)1 <3. (3.21) 
From (3.20), (3.21) and (3.2) we obtain 


i] Cy a ailie 
ba On On 
where C, is a constant dependent on M. 
We transform the left-hand side of (3.19) as follows 


f a(x)v(x)dx= f a(x) div[W(x)(x—n)] dx 


C C 
S53 eo (2) S75 81 (3.22) 


= Ja[W-(x—n)],do— J a: [W-(x—n)],do 
aa ay 
+ J «[W-(x—n)],do— J «-[W-(x—n)],do 
ie ae 
+f @[W-@—ndo-— fo [W-—n)],do 
Zha Zia 
+ J a [W-(x—n)],do— f a [W-(x—n)],do. (3.23) 
2H a 


where [W. (x — )], is the component of the vector W (x) . (x — y) normal 
to Landd> 4h. 
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Let us divide the terms of the right-hand side of (3.23) into three groups 
and estimate each group separately. 

Let us consider the first group. By virtue of Lemma 3 the function W (x) 
on the surface oY, Dive differs from +1, — 1 respectively, by no more than 


h \% 
Cs; (+) ’ 
A 
where C; are constants dependent on M. Therefore 


J «-[W-(x—n)],do- J a-[W-(x—n)],do+ [ «a [W-(x—n)],do 


Py py a 
— J «-[W-(x-n)],do 
s 
- h \“ 
> w(Qz)+ u(Qz)-—Cy (+) (3.24) 


It follows from (3.2) and the definition of X), that if the straight line 
defined by the angles ¢, 6 intersects a , then the following inequality holds 


If: (G, 8)— fo (9, MISCgA. (3.25) 
It clearly follows from (3.25) that 
p(Q)—p(Q7)—B(Qz)SCy°A (3.26) 
where Cy depends on M. Substituting (3.26) in (3.24) we obtain 


J a-[W-(x—n)],do— J « [W-(x—-n)],do 


Dy pts 
or J a-[W-(x—n)],do— J «-[W-(x—-n)],do 
zy zy 
h \° 


We consider the second group. Clearly 


J a-[W-(x—n)],do— J a*[W-(x—n)],do <a (2) ,)+0(241)- 


i 
Lha Ln 


(3.28 
We consider the function 


wo(1) StI = EOE) 


a(2')+a(2') 
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The function w@ (A) satisfies the conditions of Lemma 4, and thus for any 
t < 1 and a there will be a number & such that 


w(t) —a(t<—, k<n. (3.29) 


(<4) (3.30) 
A=t" 


Substituting (3.30) and (3.29) in (3.28) we obtain 


We set 
h = t* +1 


J a (W-(x—-n)],do- ) a: [W-(x—n)] do aCe (3.31) 


i i 
Lha xn 


Finally we consider the third group. It follows from Lemma 3 that the 
function W (x) on the surface 2/, 37 satisfies the inequality 


h \Ow c 
IW (x)1<Cyy i =C,,t", (3.32) 


L 


from which it follows that 


J a-[W:(x—n)],do— J a-[W-(x—n)],do<C,3t%". (3.33) 
zi zi, 
Substituting (3.27), (3.31), (3.32) in (3.23) and substituting (3.22) and (3.23) 
in (3.19) we obtain 


C Cae 
H(Q)S Cyt + Cott +12 + CBO te (3.34) 
In (3.34) we set 


t=— (3.35) 


where 7 satisfies the relation 
1 2n+} 1 2n-1 
gatas ME ep ; ; 
Substituting (3.35) and (3.36) in (3.34) we obtain 


C 
p(Q) ST 


from which, by virtue of (3.2), the inequality (3.3) follows. 
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§ 2. A Class of Nonlinear Integral Equations 


Let us consider the equation 
Ag=f (3.37) 
where ~, f are continuous functions of x on the interval [0, 1] and the operator 
A can be represented in the form 


1 (8) 
Ag=f f§ P(x,é,m)dndé (3.38) 


0 90 


with P (x, €, 7) continuous and non-negative. The operator (3.38) is a special 
case of the URYSOHN operator 


Ag= JK (x,¢,9(6))d¢ 


when the kernel K (x, &, 7) is a continuously differentiable function, satisfying 
the condition 


K(x, €,0)=0 


Let us formulate a uniqueness theorem for equation (3.37) for a class of 
operators A having the representation (3.38). 


Theorem. Let the operator A satisfy the following conditions. 


1. For any @;,2,Q> (x) => 9, (x) satisfying a LipscHiTz condition with 
some constant C, the operator 7 Ag, exists where A, is the FRECHET deri- 
vative of the operator A at the point g. Moreover, for any w satisfying a Lip- 
SCHITZ condition with the constant C, the function A>) A, W satisfies a Lip- 
SCHITZ condition with the constant C, depending on C. 

2. If pz (x) > @, (x) the operator AZ) A,,, can be represented as 


: 2 1 
Aj, Ag,W = J Qo, EZ g, P2 (€)-—, (é)] w(é) dé 


where the function Oo, (x, &, 7) is non-negative and continuous for 7 > 0. 
3. The integral 


11 
J J Q,(, én) dédn 
0 0 


is convergent for any g and x. 
4. There exists a continuous function of two variables w (t, h) dependent 
only on the constant C and satisfying the conditions, 
a) w(t, h) > 0, kh > 0, t+ 0; 
b) for any ¢ satisfying a LipscHiTz condition with the constant C the ine- 
quality 
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hQ, (x, ¢, h) 
W, (x, h) <w(x-€, h) 


holds where 
Axth 


W, (x, h= J JQ (x, &,h) dédy. 


Then the solution of equation (3.37) is unique within the class of functions 
satisfying a LIPSCHITZ condition. 

It should be noted that this type of conditions is satisfied by the integral 
operators arising in the inverse problems of potential theory. 

The classical example of an inverse problem is the one arising in the theory 
of NEWTONIAN potential formulated by Novikov. Solving the planar inverse 
problem of Novikov is equivalent to solving equation (3.37) with 


P(x, €,h)=In [r? +h? —2rnycos22(x—£)-n], O0<o<r. 


If ~, (x), @2 (x), and w (x) are continuous functions, @, (x) > @, (x), then 
in this case the function A; A,,W is a density of a single layer potential with 
weight @, (x) distributed on the curve defined by the polar equation 


Q= P(x) (2) 
(2 xx is the polar angle) and satisfying the following condition, the logarithmic 
potential of a single layer with density @, (x) w (x) on the curve 

e=9,(x) (W,) 


coincides with the single layer potential with density Ae A,,w outside the 
domain bounded by J°,. 

It is easy to show that in this case the operator Aj. A 4, satisfies the condi- 
tions 2, 3, 4 of the theorem. 

Let us introduce some notation. Let g; be continuous functions and let 


u (x, P1> 2) =max p(x) ’ 
j=1,2 
v(x, Py, @2)=min ~, (x), 


W(X, 91,92,93)= min w(x,9,9)= max v(x, ¢,). 
i, j=1,2,3 i,j=1,2,3 


Leta; (j = 1, ...) be sets of points of the interval [0, 1]. We denote by 0 (1, %) 
the distance between the sets «,, «5. 
We denote by « (@,, @,) the set of values of x for which 


P1(X)> 2 (x) 
and by & (~,,@,) the set of values of x for which 


P1 (x)=@2(x). 
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Let A; j = 1, ...) be some real numbers. We denote by « (~1,@>, 4) the set 
defined by the relations 


: XEO(P1, 92,4) 
| 


XEA(D1, Pa), oLx,%(91,92)]=4 
and by & (@;,@5, 4) the set 


&(P1, Pr, AV=24(P1, 2) OE (G1, Y2)Oa(P1, G2, A) 
while by &(@;,@>, 4,, A.) the set 


&(Q1, G2, 44,42) =A (Py, P2,4,)OX (Gy, Pz, 42). 


The proof of the theorem is based on the following lemma. 


Lemma. Let A be an operator satisfying the conditions of the theorem and 
Iet p, andq@, be functions satisfying a LipscHiTz condition with the constant C, 
and such that the sets « (~,,@2) are not empty. Then for any sufficiently small 
positive A and / there exist functions @,, @,, W possessing the following prop- 
erties, 


1. U(X, Gy, Pa) 2G; (X)Z(X, 1,92), JHl,2; 
2. Ag,—Ag,=AG,—AG, + AW, U=U(x, Py, G2); 
3. W(x)20, xEa(P2, 9): 
W(x)S0, x€a(P;,92); 
4. e[a(P1, 92), “(2,G,)|>A; 
max |$2(x)— (I=; 


6. the functions W, G; (7 = 1, 2) satisfy a Lipscuitz condition with the 
constant C, depending on C. 


We prove the above lemma. Let @,,@, be functions satisfying a LipscHITz 
condition with the constant C, let 


Ag,—Ag,=f 


and let 2, 4 be some sufficiently small constants. 
We denote by y,,, the function defined by the relations 


1 x 
a) Yma=—— (2 —1)s XE4(P2, P1,4);3 


1 7 
b) Yma7?— (P1 — P2)s xEE(P1, P24); 


Yma= 0, XEL(Py, P2, 2A Ox (G2, G1, 2A); 
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c) Ym is a continuous function on the interval [0, 1) and it changes linearly 
on the set &(@,,@2, 2 4, A). 
We consider the sequences of a function Q mrasP2mna WV inka (A = 1, -- Dy 
P jmoa = Qj ! (j=1,2) 
D jmk+12=P jmkat Ojmka 


W mod =0, Win, k+1,4— Winka + Omka ? 
where 


a) Stina = mas Pamka— Pimea)s xEa(2, 94), 
5 2mka= mas Pimka— Paka)» XEX(Py, Pr); 

b) Sima =W [Arad SiO meas 0; Pamka—Pimeal> X€x(P1, 2), 
S2mka= W [Abinics 4 omicaD mia 5 03 Pima Prmeal, XEa(P2,91); 
Umea (X)= U(X, Prmkas Pamka) » 


Umka(X) = V(X, Pi mkas Pramna) » 


Sina (X) = 5 a mua (%); XEX(P2, P41); 
Sinz (X) = Samta (X), XE4(P1, 92); 
C)  Ojma=9, XEX(Py, P2)3 
4) Onxa=A ae vaca mk —Oimnas XEX(P1, 2), 
Onka= Avra dom waOrnk —Oamka» XEA(P2, 91), 
Sint = Aurmna4 oe XEL(Py, Pr). 


One can easily see that 


C,k 1 
[A@2mka— AP meat AumcaY mka— sis, (4) (3.39) 


where C;, is a constant depending on C and y (f) is the modulus of continuity 
of the function P (x, &, 7). 
Indeed by virtue of the definition of the FRECHET derivative 


AQ jmk+1a— AQ jmka— A oo jinx imkall| =0 (15 jmzalD , (3.40) 


From the representation (3.38) of the operator A it follows that the right- 
hand side of (3.40) satisfies the inequality 


0(I5 pals (4). 


The inequality (3.39) now follows by virtue of the definitions of the functions 
Pj mka> Oj mka W mae 
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Next, from the conditions imposed on the operator Ay Ag, it follows that 


for k < m the functions Qj mya, Wma Satisfy a Lipscuirz condition with a 
constant C3 dependent only on C and that fork > m 


P jm, k+ 1,a=P jmca=Pjma> 
~ (3.41) 
Wn, k+1,2=Wmka=Wma- 


The functions @;,,, in (3.41) evidently satisfy the conditions 


O[« (Pima Prma)> ~ (Pama Pima lZ; 


i 3.42 
0(Q1, 92) SP jma SU (G1; Yr) Ae) 

and for any sufficiently small , A 
Max |P2ma—Pimal > 2h. (3.43) 


x 
We consider now the sequences of functions G, mas Bomar W ma (m = 1,...). 
By virtue of the fact that all the elements of sequences are bounded and satisfy 
a LipscHitTz condition with the constant C, independent of m, one may choose 
from the sequences uniformly convergent subsequences 


Pimgas Pamga Winga (q=1,...). 
By virtue of (3.42) and (3.43) the functions 
@,;= lim P jmaao (j=1, 2) 
i 
17 = lim Wad 
q7o 


satisfy the following relations 


e [a (G1, G2), «G2, OiJ2A, 

AG,—AGz +AW=f, 

v(P1, P2) SG; SU (G1; 2), (3.44) 
max |G.—6,|>h 


x 
and a Lipscuirz condition with constant C4. 
We denote by y,, the function defined by the relations 


Pee, 2 2 oe 
Ym=— (G2 — 1), xEa(G2,G1), 
fine | ad os aa 
Pm =— (G1 — G2)> x€a (G1, 2) 


Vm =O x€&(G1, G2) 
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and consider the sequence of functions @ mx, P2mme YW me (kK = 1,.-.) 


QP jmo= Pj (j=1,2), 
P jm, k+ 1 =P jmk tO jmk > 
WVno=W 5 


Win, e+ = Wink + Omk » 


where 0; mx 5 my are functions defined by relations similar to relations a), b), c) 
defining the functions 6; mkas Oma 
In complete analogy with (3.39) we obtain 


Ck {1 
|AQ@omx~— AP imk + Aum mk — JF | <S,(1) ’ (3.45) 
max \P2mp— Pimpl >h ? 
MAX |P2mp+1—Pimptil Sh. 
By virtue of the definition of the sequences (;,,, we can obviously find a 
number p, 0 < p < m, satisfying (3.46). 
We denote by G;,,, W,, the functions D,m = Pjmp» Wm = W mp and consider 
the sequences 


(3.46) 


Pim Wm (m=1,...). 

In analogy with the sequences 2; mj, Wing (m= 1,...) we may select uni- 
formly convergent sub-sequences Pimg? 17 mg (q = 1,...) from the sequences 
Dim W m- By virtue of the above, the functions 

P;=liM Pim, = (F= 1,2), 


q7 a 
v=lim ¥,,, 
q7o 
satisfy the relations 1)—6) and the lemma is proved. 

We now prove the theorem. We suppose the solution of (3.37) to be non- 
unique within the class of LipscHITz continuous functions, i.e., there exist 
distinct functions @,,@, satisfying a LipscHiTz condition with the constant C 
and 


Ag,—Ag,=0. (3.47) 
It follows from the Lemma that for any sufficiently small A and A there exist 


functions @,, @., W satisfying the conditions 1)—6) of the Lemma and there- 
fore, according to (3.47), 


AG,—AG,+A,=0. 


Applying the operator AS to (3.48) and making use of the conditions 2) 
and 3) of the theorem, we obtain 
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1 
J J OaLx,¢,n] dn-sign[ 2 (€)—G, ()ydé+P(x)=0. (3.49) 
We denote by X a value of the independent variable x such that 


192 (X)—@, (*)|=max |P, (x)—G, (x). 
Because of the properties of @; . 
192 (X)—@, (X)l=h. (3.50) 
Let, for the sake of definiteness, 
P2(X)>@, (X), 
1.€., 
P2(X)—G,(X)Hh. (3.51) 
We put x = ¥ in (3.49) and divide the integral in the left-hand side of (3.49) 
into three parts 


X+hy a (6)—v (4) ; 
Q(X, 6,1) dy sign [G2 (2)—G, (E)} dé 
xX—hy 0 


¥—hy @(S)-0 (8) 


+ fj J 0;(X, €,h) dy-sign[@,—@,J} dé (3.52) 


0) 
1 ae) 
+ J § Q,(%,¢,h) dn-sign[@,—9,]dé+P(x)=0, 
x (0) 


Ay 
where 


h,= Z 
a oe 

We denote the integrals on the left side of (3.52) by J,, J, and J; respectively 

and estimate them separately. Let us consider the first integral. From (3.51) and 

condition 6) of the lemma it follows that 


h 
2O-BOZF — s-FlEye, (3.53) 


from which 
1, >W,(3,h), 


where W, (x, h) is the function defined in condition 4) of the theorem. 
We consider now the second and third integrals. From condition 4) of the 
lemma we have, 


X-AG(S)—v (5) 


1,2 i) ) 0;(x; €, n)dy-sign[@,—-@,] dé, (3.54) 
0 @) 
1 a(&)—v (6) ; _ 2 ; 
=f  f§ Q,(%,é.n)dn-sign[9,—9,] de, 
X+A 0 
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while from (3.54), condition 4) of the theorem and condition 5) of the lemma 
we obtain 
I,4+1,>—-—o(/A,h)-W,(x,h). (3.55) 


From condition 3) of the lemma and from (3.51) it follows that the last term 
in the left-hand side of (3.52) is non-negative, i.e., 


W(x)>0. 


From (3.54), (3.55) and (3.56) it follows that for sufficiently small A there holds 
the inequality 
I, +1,4+1,+(x)>0 


contradicting the equality (3.52). Thus the theorem is proved. 


§ 3. Inverse Problems for Some Non-Newtonian Potentials 


Inverse problems for NEWTONIAN potentials have unique solutions only 
under very essential restrictions on the nature of the density distribution. An 
example constructed by P. S. Novikov indicates that the solution of the inverse 
problem of NEWTONIAN potential may be non-unique even in the class of 
densities which are piecewise constant with a given constant and are contained 
in simply connected domains. An analogous situation seems to take place 
also for potentials corresponding to the fundamental solutions of other elliptic 
equations. 

As will be shown in this section, for inverse problems with some Non-NEw- 
TONIAN potentials not connected with differential equations, uniqueness holds 
for much milder restrictions on the nature of the density; its continuity and 
finiteness will be quite sufficient. 

Let x, € be n-dimensional vectors with components (x, ..., Xp), (&) «5 En)s 
let D, D, be non-intersecting simply connected bounded domains of n-dimen- 
sional] space and let 


: OoaG é) 


“o(s)de, (3.56) 


where 
r(x, g)=|x—¢]. 


Given the function u(x) in D,, it is required to reconstruct the potential 
density oe in D. 


Theorem. For p > 1 + 2 the solution of the problem is unique in the class 
of continuous @ (x), i.e., the continuous function 0 (x) is uniquely determined 
by the values of the function u (x) in the domain D,. 
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For our proof we shall examine the function 


1 
v(x) = | —— o(é)dé S57 
Sao (é) (3.57) 
where 
: ye Cet In Le 
R(x, 9, 0) = Jr? (x, 2) +1? 
and y is a vector with components (),,..., ¥m41)) MI= p — uw — 3. 


The function v (x, y) will clearly be a harmonic function of the variables 
x, y with cylindrical symmetry. From the results of § 6, Chapter II it follows 
that the function v (x, y) is uniquely determined on a set Q of x, y space: 


x, yEQ, if XED or y+0; 
by its values on the set 
xeD,, y=0. 


The function v(x, y) is a NEWTONIAN potential of a simple layer with 
density 9 (x) distributed on the n-dimensional manifold x éD, y = 0 in the 
(1 + m+ 1)-dimensional x, y space. It is clear that 9 (x) is uniquely deter- 
mined from v (x, y). Hence the theorem is proved. 

We note that the method used for the proof of the uniqueness theorcm 
permits one to obtain estimates characterizing thc stability of the solution. 


§ 4. An Inverse Problem for the Wave Equation 


In this section we shall prove a theorem giving the uniqueness for the 
solution of an inverse problem for the wave equation. 
We consider the equation 


2 
gee a (3.58) 


where uw is a function of the three variables x, y, ¢ and # is a function of the 
variables x, y. 

We consider the following problem: the domain Dg is given in the x, y 
plane. The function n(x, y) > 0 is continuous and identically equal to one 
outside of Dj. Moreover, in some domain D, with D, MO Dy empty, there is 
given a family G of solutions of (3.58) for all ¢ > 0. It is required to determine 
the function 7 (x, y) inside Do. 

We note that a similar problem in the case of one variable was considered 
in the paper [59]; quite close to the problem being considered there, is the well- 
known STURM-LIOUVILLE inverse problem. 

We formulate a uniqueness theorem for the above problem for a certain 
family G. 
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Theorem. Let D, and D, be bounded and simply connected domains, and 
D, a bounded, simply connected domain, which does not intersect with Dy or 
D,. Let G be a set of solutions of (3.58) satisfying the following initial condi- 
tions 


u(x, y,0)=0; 
0 
£ u(x, ¥,0)=5 (059 Yo) (3.59) 


where Q (Xp, ¥o) is any point of D,. 

Then the solution of the inverse problem is unique, i.e., the function 
n (x, y) is uniquely determined within Do. 

We give a brief sketch of the proof. We denote the solution of the CAUCHY 
problem (3.59) for equation (3.58) by u (x, y, Xo, Yo, ¢) and consider the function 


v(x, Y,Xo, Vo A= f U(X, Y,X9, Vo, t)cosAtdt. 
One can easily see that ; 
Av= —5(x—x9, y— Yo) —A?nv. (3.60) 
The function v is the fundamental solution of the HELMHOLTZ equation 
Av=—J*n?v 


with the singularity at the point Q (%, yo). It is well-known that a fundamental 
solution for an elliptic equation with analytic coefficients is an analytic function 
of both the independent variables x and y and the coordinates of the singularity 
Xo, Yo everywhere away from the singularity. 

According to the conditions of the theorem the function v is given in the 
domain D,, of the 4-dimensional space R (x, y, Xo, Yo), the direct product of 
the domains D,, D, of the spaces P (x, y) and Q (x, yo). The function 7 (x, y) 
is identically equal to one outside of Dy. Hence, by virtue of the uniqueness of 
the analytical continuation, the function v may be considered to be given 
everywhere outside of the domain Dog = Dy X Do. 

Let Dg be a bounded domain containing the domains Dp, D, and D, and let 


P(x +(y—n)’, ro=(%o—8) +o)’. 
From (3.60) it follows that for Pe Dg 


1 
v(x, Y>Xos Vos A=5_ In [(x—x9)?+(y— yo)" ] 


1 
—5— An? (x,y) i) 0(6,9, Xo, Yo' A) nr dé dy 
rh Do 


+5(X, y, X95 Yo, €)3 (3.61) 
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where J’3 is the boundary of Ds. 
We denote by 2, (x, y, Xo, ¥p) the function 


D cA 0°v(x, YsXo. Vo, 9) O7v(X, ¥,Xq, Yo: 0) 
= + : 
OA OX 0X5 Oy: Oyo 


[t can be easily shown that the function v,, for RE Ds3 = Ds % Dy, equals 


2 J me ee NOD ag -dnt 3, (3.62) 


where 


ee sje © 075 i 075 
I X, ¥,%o, Jo ah Ox OX, Op-Oyo : 
By virtue of (3.62) the functions v,, 5, for RE Dgg and RE Dyy are analytic 
functions of the variables x, y, xo, Vo. 
By setting 
X=NXo, Y=Yo; P(x, y)eD, 
in (3.62) we obtain 


7 1 ‘ 
D2 =(v, (x, VX, y)-v, (x, yx, y= ) n (¢, n) a abdi 3 (3.63) 
Do 
The assertion of the theorem follows immediately from the results of § 3 
and (3.63). 


We note that the theorem generalizes to the case of the wave equation in a 
space of an arbitrary number of dimensions as well as to the heat equation and 
some equations of hyperbolic and elliptic type of higher orders. 


§ 5. On a Class of Inverse Problems for Differential Equations 


In this section a statement more general than that of the previous section 
is considered. A class of inverse problems is reduced to linear integral equations 
of the first kind which are then investigated by methods differing from those 
of the previous section. 

We consider the equation 


) ) 
Py (2) u(x, y)= Ps (2) P, (2) u(x, y), (3.64) 


where x, y are vectors with components (x,,...; Xp), 4... Yq) and P; are 
polynomials with coefficients continuously depending on x. 
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In regard to P,, we assume that the equation 


P, (2 )oco=0 


has a fundamental solution G (x, x°) defined in the entire x, x, space*. 

We pose the inverse problem for equation (3.64) which will be considered 
in this section. 

Let Do, D, be some bounded domains of x space, where the intersection 
D, © Dg is empty. The coefficients of the polynomials P,, P, are given for the 
entire x space and the coefficients of the polynomial P3 are given everywhere 
outside of Do. In addition, let a family of solutions u (x, y, €) of (3.64) depending 
on the parameter é (£,, ..., €,), be given in the domain D,. 

The functions w (x, y, §) satisfy the conditions 


r 


— © u(x, 0, C= Ser ++ km (% S)> 


Oytt... Oym 
O"'u (x, 0, &)| JOG (x, x° 
my eeken chic A Zee ee (3.65) 
a7 ow Sexo Ox} Save Ox," Ox}? bai Ox} 
Ou (x, y, ) 
eee <C,+lyl™ 
Oy; aed On" 
where the indices ,, ..., 4, rum over the values which are less than the maxi- 


mum values of the corresponding indices in the polynomial 


0 
*3 (=) 


and i; + j,,...,/, + J, run over a similar set with respect to the polynomial 


0 


pee is a sphere of radius 0 with the center at the point x°, f(x, &) are given 
finite functions, and C,, C, are some constants. 

It is required to determine the coefficients of the polynomial Pg, in the 
domain Dp. 

In the case when 7 = #1 = 1 the above inverse problem is equivalent to 
the known SturM-LIOUVILLE inverse problem for ordinary differential equa- 


tions. 


* The existence of a fundamental solution has been proved, say, for elliptic equations 
of the second order and also for ordinary regular equations with constant co- 
efficients. 
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We denote by v (x, A, &) the function 


v(x, A,é)= J u(x, y, €exp {—(A, y)} dy 


(A is a vector with components A,, ..., 2,,)- 


By virtue of (3.64) and (3.65) the function v (x, A, &) satisfies the following 
differential equation 


3 
(2 )omry(—2y-Pa(2 Jory. (3.66) 


J Jj 


where 
0.6 (OO 
x,A4,€é)=O| A, — ]}-u(x,0,€ 
W(x, 4,8) o( i Be, 2) (x,0,€) 
is a function which appears as a result of integrating by parts the expression 


for 
0 
P,{ — ]}- 
: (=) : 
and Q is a polynomial. 


We multiply both sides of (3.66) by the function G (x, x°) and integrate 
over the whole x space. From inequalities (3.65) and (3.66) it follows that 
the integrals will converge, and by integrating the corresponding expressions 
by parts, the boundary terms at infinity will be equal to zero. 

Thus, we find that the function v (x°, 2, &) satisfies the integral relation 


v(x°,2,2)= J G(x, x°) P34) Pa (2) vdx-+ J G(x,x°) pdx. (3.67) 


From (3.67) it follows that the function v (x, 4, €) is an analytic function 
of A in some neighborhood of the origin. We consider 


r 


3 
v reek KC =——_—_— v(x, 0, ¢), 
kao 0 8) OAK. Ofte 


r 


0 
Wr, rec =e pie O09): 


pts Bake Ode 
By virtue of (3.67) the functions vy, .-- 1 Whey +++ kyo Pky «+ km Satisfy the 
following system of relations 


Vikey pe s)= { G(x, ae es a (2)x Vivtnis coh dx 


+ § 60% x) Vag a OO AX, (3.68) 
Vp 9 (X89, = J G(x, x") Wo. 0 (x, 6) dx. 


P, P;(0). 


BeatriceGloria_personal library 
68 Inverse Problems for Differential Equations 


Relations (3.68) may be regarded as a recursive system of linear integral 
equations of the first kind in Dg for the coefficients P;,... ; of the polynomial 
P,. In fact it follows from (3.68) that the function ug... aes &) may be re- 
garded as given in the whole space. 

We consider now the function 


re) 
Vo..1 ees €)= J Gi, x°) Po 1... 022 a— ] 20... 0%, 6) dx 
—=——_ Do —_— e) 
q q 


J 


re) 
a J G(x, x°) Po “nil 0Pa( 2) 06 wold 


Do = i 


Bs A) G(x, x") Wo ...1...0(%, ) dx. (3.69) 

a 
If x° € D, then in the equality (3.69) all the terms may be regarded as given 
with the exception of the first one in the right-hand side, and therefore (3.69) 
may be regarded as a linear integral equation of the first kind for the function 


Prong hig): xEDpo 
—— 


q 
with the kernel 


K(x, oe O=G(x, x) P, (2) Vo 0%, g), xeD, ° 


If the above equation has a unique solution for any g then the equations 

defining the function 
Po 1 Oe ing x) 
—_ 
q qi 

is found in a similar way, and so on. 

We will now carry out a detailed consideration of the system (3.68) in the 
case in which the original differential equation is of the form 


B J 
Aru (x, y=), a; (5m (x, y) (3.70) 


(y is a scalar, «, B are integers.) 
The system (3.68) in the case under consideration will take the form 


Op (x°, C= i} G(x, x) Ug (x, €) a, (x) dx, 


V9 (x°, S)= J G(x, x°) Wo (x, E) dx, 
k-1 


6, (x°, = —- d J G(x, x°) og - 5 (x, 8) a(x) dx +0, (x°, 6) (3.71) 


aa G(x, x°) a(x) 09 (x, E)dx— § G(x, x°)-W, (x, &) dx, 
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yf cee oad tea 9 mee for n even, 2a—n>0; 
G(x,x°)= 


id boa 


Yan |X—X in the remaining case. 


From (3.71) it follows that for uniqueness of the solution of the inverse 
problem under consideration it is sufficient that there be a unique solution of 
the integral equation of the first kind with respect to the function a (x) 


J G(x, x°)+v9 (x, &) a(x) dx =i(x°, €). (3.72) 


We multiply both sides of (3.72) by an arbitrary function p (x°) and integrate 
over the domain D,. 


Jo(x°, 8) oe) dx°= J W(x) v9 (x, Ga(x)ay, 


Di 


W(x)= J G(x, x°) a(x) dx?®. (3.73) 
Dy 


The function W (x) in (3.73) is a potential with density 9 (x9) distributed 
in the domain D, whereas the function vg (x, §) is a potential with the density 
Wo (x, §) distributed in the domain D,. Hence, for uniqueness of a solution 
of (3.72) it is sufficient that the linear hull of products of these potentials should 
be dense in the set of functions defined in Dp. 

Let, for the moment, the set of functions Wo (x, ) be such that the set of 
potentials vy (x, €) is dense in the set of all harmonic functions regular in some 
extension of the domain Dy — Do, *. 

From potential theory it is known that the potentials W(x) lie densely in 
the set of all harmonic functions regular in Dp,. It can be easily seen that in 
this case the linear hull of defined products of potentials is dense in the set of 
continuous functions defined in Dy. In fact, we consider the harmonic poly- 
nomials of order v which are normalized and homogeneous with respect to 
the coordinates of the points x — x! where x! is an arbitrary point of Dp. 
It is well-known that the sum of squares of all these polynomials is equal 
to |x — x1|* whereas the function {1 — |x — x*|?/R?} approximates 6 (x — x1) 
in the sphere |x — x!| < R. Thus, we have proved the following uniqueness 
theorem for the inverse problem for equation (3.70). 


Theorem. If a family of solutions u (x, y, £) of (3.70) is such that the set 
of potentials vy (x, €) is dense in the set of all harmonic functions regular in 
the domain Dg, (some extension of the domain D,), then the solution of the 
inverse problem for (3.70) is unique in the class of continuous functions a, (x) 


(2s Peery oh 
* This condition is fulfilled if, for instance 
Po (x, &) = u (x, 0, &) =4 5d (x— §) 


where the parameter & runs over the domain D,. 
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